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Abstract 



We review SO (10) grand unified theories (GUTs) in four and five dimensions (4D and 5D). 
The renormalizable minimal SO(10) SUSY GUT is the central theme of this review. It is 
very predictive and makes it possible to construct all mass matrices including those of the 
Dirac and heavy right-handed Majorana neutrinos. So it is not only able to reproduce all the 
low energy data, except for too larger # 13 in the lepton mixing angles (which can be evaded 
without spoiling the basic ingredients), but also predicts almost all new physics beyond 
the standard model (SM) like neutrinoless double beta decay, the electric and magnetic 
dipole moments of the quark-lepton, lepton flavour violation, leptogenesis etc. To be very 
predictive, on the other hand, implies that predictions are unambiguous and they are always 
exposed to severe compatibilities with observations as well as to a conceptual consistency 
check. The explicit construction of the Higgs superpotential and the explicit display of 
a symmetry breaking pattern from GUT to the SM show that the naive desert from the 
SM to GUT in the minimal supersymmetric standard model (MSSM) is a too simplified 
concept and we have many definitely determined intermediate energy scales in general. This 
situation destroys the naive gauge coupling unification in the MSSM scheme. Also the 
precise measurements of neutrino oscillation data have revealed several small but manifest 
mismatches with our predictions. Also there are arguments that it is impossible to construct 
a GUT theory in 4D with a finite number of multiplets that leads to the MSSM with a 
residual R symmetry. Of course, there are several loopholes which solve these problems in 
4D. However, if we try to solve all these pathologies comprehensively, it is very attractive 
for us to go into extra dimensions. Extra dimension may be either warped or flat. The fifth 
dimension, for simplicity, is compactified on the S 1 /( y Z 2 ) (warped) or on the S 1 /(Z 2 x Z' 2 ) 
(flat) orbifold with two inequivalent branes at the orbifold fixed points. In the former warped 
case, intermediate energy scales are translated with the positions of Higgs fields in the bulk 
and the fundamental scheme of the MSSM is recovered. On the other hand, in the latter 
flat scenario, all matter and Higgs multiplets reside on the Pati-Salam (PS) brane where the 
PS symmetry is manifest. There the original renormalizablity in Yukawa coupling is broken 
but its essential structures of mass matrices in minimal SO(10) GUT in 4D is promoted 
to the PS invariant action in 4D. In the gaugino mediation mechanism, the SO(10) gauge 
multiplet is transmitted to the PS brane through the gaugino mediation with bulk gauge 
multiplet. Further breaking of the PS gauge group to the SM group is realized by VEVs 
of the Higgs multiplets (4,1,2) © (4,1,2). We show that this model not only cures all 
pathologies in SO(10) GUT in 4D but also provides the consistent inflation scenario and dark 
matter candidate and leptogenesis. The gauge coupling unification is successfully realized 
at Mqut — 4.6 x 10 17 GeV after incorporating the threshold corrections of the Kaluza- 
Klein modes, with the compactification scale (assumed to be the same as the PS symmetry 
breaking scale) M c = v-p$ = 1.2 x 10 16 GeV. This orbifold GUT model can naturally leads 
us to the smooth hybrid inflation, which turns out to be consistent with the WMAP 5-years 
and 7-years data with the predicted Mqut and vps- However, this simple model suffers from 



the stau LSP (the lightest SUSY particle) problem. Neutralino LSP can be realized when 
the compactification scale of the fifth dimension is higher than the PS symmetry breaking 
scale, keeping the gauge coupling unification. We reanalyze all new physics beyond the SM 
in 5D, the leptogenesis, LFV etc. So our predictions range over all particle physics. The 
recent discoveries of a Higgs-like object and null search of LHC give very clear and important 
inmpacts on the above mentioned scenario and we discuss to where they drive our theory. 
Finally we add some comments on the impacts of the discovery of a Higgs-like particle by 
the Large Hadronic Collider at CERN (LHC). 
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Chapter 1 
Introduction 



SUSY GUT is the most promising model beyond the SM |IJ. The SM is a very powerful 
theory but it has validity limits like the other great theories. There are discrepancies with 
observations like neutrino mass as well as the other non affirmative ones like muon g-2 @,an 
anomalous like-sign dimuon charge asymmetry ||, and h — > 2j. The CKM phase only is 
insufficient for baryogenesis ||. The SM does not predict dark matter (DM) candidates. 
These facts strongly suggest a more comprehensive theory which improve the deficits of the 
SM. 

Above the above mentioned observational deficiencies, the SM has conceptual problems 
which do not allow it to be complete. Indeed, the SM+z/ has too many parameters, 19 + 
9. It does not explain quark-lepton mass hierarchy, mixing angles, CP phases, Higgs mass 
stability against quantum corrections, three different gauge couplings and their unification. 
They are all independent given or unknown objects and have no mutual relationships in the 
SM framework. 

There are the motivations for us to consider the theory beyond the SM from the bottom up 
approach. The top-down approach from string theory is also necessary and complimentary 
to the bottom-up approaches U |j| . 

Even if we restrict ourselves in the Higgs hierarchy problem, there are many approaches 
like little Higgs 0, composite particle models of technicolor || etc. However, it seems to be 
only SUSY GUT which may solve all these problems comprehensively. 
MSSM explains hierarchy problems and gauge coupling unification. However, it gives rather 
few definite predictions beyond the SM since it gives no information on mass matrices of 
quarks-leptons and on the physical structures at the unification scale. This drives us to 
SUSY GUT. Even if we accept a scheme of SUSY GUT, there are so many options. What is 
the gauge group, SU(5), SO(10), E 6 , E 7 , E$ or their products ? So we need the additional 
criteria to select the gauge symmetry. An anomaly free condition may be a good candidate 
for it since chiral symmetry must be preserved under quantum correction. SO(10) is the 
smallest group which is free from anomalies in a single multiplet. Such an anomaly free 
condition is meaningful if the theory is renormalizable. Even if we fixed a gauge group, we 
need other criteria to make a definitive model. Renormalizability will be one such criteria. 
Before proceeding to the detailed discussion, let us consider the structure of the existing 
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theories. For SU(3) x U(l) theory at fi < (9(10 2 GeV), it has the following form 

L = L ren + ^- (1.1) 
in general. Here the first L ren denotes renormalizable Lagrangian and the second term an 

1/2 

unrenormalizable effective Lagrangian. L\ implies the fermi coupling oc J^J^ and A x = 
O(10 2 GeV) in the effective weak interaction before the SM. For SU(3) x SU(2) x 17(1) SM+z/, 
this L\ becomes the renormalizable g2J fl W fl term but a new effective term appears in the 
(type I) seesaw mechanism ||, 

L = L' ren + Yj^-Y v {LHf ee L' ren + (1.2) 
M R A 2 

Here L' ren D L ren + jj^ and A 2 (= O(l0 13 GeV)) ^> A^ 2 . Thus the theory is expressed as 
the sum of a renormalizable theory plus effective action with cut off, and the renormalizable 
Lagrangian becomes more involved as the energy scale goes up. When the energy scale n 
goes up and // > Mr, this Lagrangian is expected to be renormalizable. This is the rough 
sketch of GUT schemes. This scenario, of course, is not a unique one. For the other types 
of seesaw mechanism [|T0[| , |TT[| , and |12| (Type II, III, and Inverse seesaw, respectively) give 



the corresponding phase transitions at each energy scale. However, even in these cases, GUT 
scheme that the renormalizable part becomes enlarged as energy scale goes up seems to be 
universally valid. Especially, minimal SO(10) GUT includes Mr in the Yukawa coupling 
with the other quark- lepton for the case of 126 Higgs field, and this field is suitable for 
the Yukawa coupling. Since after GUT all local field contents becomes massless, it is very 
natural to consider renormalizability as a guiding principle for constructing GUT theory. 

Quarks and leptons have analogous structures. They have both three families of left- 
handed doublets and three right-handed singlets. However, they have quite different mass 
hierarchies and mixing angles. The successful gauge coupling unification of the minimal su- 
persymmetric standard model (MSSM) strongly supports the emergence of a supersymmetric 
(SUSY) GUT around M GUT ~ 2 x 10 16 GeV. 

However, the MSSM does not specify the gauge group at M GUT , and therefore the sym- 
metry breaking pattern to the SM gauge group. 

Group theoretical properties give strong constraints on seemingly quite different objects 
of quark- lepton. The minimal SU(5) model JL3| is ruled out from the wrong mass matrix 
relation and fast proton decay [Q, and, among possible gauge groups, SO(10) GUT seems 
to be most promising. SO(10) |]T5| is the smallest simple gauge group under which the entire 
SM matter contents of each generation are unified into a single anomaly-free irreducible 
representation: 16 representation. This 16 representation includes the right-handed neutrino 
and no other exotic matter particles. Among several models based on the gauge group 
SO(10), the renormalizable minimal SUSY SO(10) model (minimal SO(10) GUT) has been 
paid a particular attention, where two Higgs multiplets {10 © 126} are utilized for the 
Yukawa couplings with matter 16j {i = generation) |16[ Jl7| |H| . 



Before discussing the details of this model, let us explain the general structure of SUSY 
GUT. It consists of three ingredients: 
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(A) To build the Yukawa coupling based on gauge symmetry. 

(B) To construct the gauge invariant Higgs superpotential and exhibit the gauge symmetry 
breaking pattern from GUT to the SM. 

(C) To show the supersymmetry breaking mechanism and the initial conditions of soft SUSY 
breaking terms at GUT . 

The SM discusses only (A). The MSSM and its variations like the NMSSM and nMSSM etc. 
deal with (A) and (C). So far minimal SO(10) GUT can uncover (B) since it can fix the 
gauge invariant Higgs superpotential. 

Let us proceed to the detailed arguments of the model. Concerning (A), a remarkable 
feature of the minimal SO(10) GUT model is its high predictivity on neutrino oscillations as 
well as reproducing charged fermion masses and mixing angles. In this review we consider 
the renormalizability as one of the principles. If we relax the renormalizability, different 
SO(10) models are also possible [|19| [EjJ |2l| etc., which will be discussed at subsection 
|2.9.3| breifly. Minimal SO(10) GUT includes the additional CP phases which solves the too 



smaller baryogenesis problem for the Cabibbo-Kobayashi-Maskawa (CKM) CP phase \T2 



However, after KamLAND data |23[ had been released, it entered to the stage of precision 
measurements, and some mismatches with our predictions were revealed in #13 in the Maki- 
Nakagawa-Sakata (MNS) matrix and neutrino mass square ratios. Recently its value has 
been measured precisely 



sin 2 2fl 13 = 0.092 ± 0M6(stat) ± 0.005(syst). (1.3) 

Many authors performed data-fitting analysis to match up these new data. Some models 
have seeked the solution of mismatches in the threshold corrections, though no one has yet 
performed these tremendous calculations. The SO(10) model has, except for the SUSY part- 
ner, no exotic matter but has many guge and scalar particles other than the SM. Such heavy 
particles or would-be Nambu-Goldstone bosons are studied in (B). The Higgs superpoten- 
tial in the minimal SO(10) GUT has been constructed and a detailed analysis of symmetry 
breaking patterns have been extensively studied by us [ER EB, E7|. This construction gives 



the vacuum expectation values (VEVs) at intermediate energy scales explicitly, which gives 
rise to a trouble in gauge coupling unification as well as the necessary scales of the seesaw 
mechanism [[| and fast proton decay. 

Probably, we have no desert between electroweak scale and GUT which MSSM has been 
assumed so far for simplicity. At each stage of the intermediate energy scales, there enters 
new massless particle into the renormalization group equation, changing the naive behaviour 
of gauge couplings and spoiling the unification. This mismatch of gauge couplings has been 
explicitly shown in Ref. ||28|| , where the couplings are not unified any more and even the 



SU(2) gauge coupling blows up far below the GUT scale. 

Even if this problem was solved, the renormalizable minimal SO(10) GUT potentially 
suffers from the problem that the gauge coupling blows up above GUT and below the Planck 
scale. The minimal SO(10) GUT also predicts too faster proton decay p5| |29| . 



So far we have discussed mainly on the observational conflicts. In order that a realistic 
GUT model works well, it must not involve the internal inconsistency concerned with (C). 
Supersymmetry and gauge symmetry are related via U(1)r symmetry. The no-go theorem 
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on SUSY breaking |30| [[yj says that consistent spontaneous SUSY breaking for gauge group 
equal and higher than SU(5) is naively impossible in 4D. 

Thus, to solve these observational and conceptual problems comprehensively, we consider 
an orbifold GUT |32| preserving the merit of the minimal SO(10) GUT. 

One of the demerits to go beyond 4D is to break the renormalizability of the preceding 
theory. However, the extra dimension opens up only in the neighbourhood of the GUT scale 
and the deviation from the renormalizable theory give only small corrections. On the other 
hand, there appear many merits by considering an extra dimension. First we can solve the 
problems mentioned above. Moreover, the problems involved in the minimal SO(10), fast 
proton decay and coupling blow up, can also be evaded. We have a new geometrical SUSY 
breaking mechanism. 

This paper is organized as follows. 

In Chapter 2 we review the minimal SO(10) GUT in 4D. In the first two sections we 
overview the historical background before we proposed SO(10) GUT model. In section 2.3 
we give a setup of the minimal SO(10) GUT model and its data fitting is given in section 2.4. 
Section 2.5 is devoted to leptogenesis and electric dipole moments (EDM). The latter is an 
important signal of new physics beyond the SM (BSM) since minimal SO(10) model has the 
definite additional CP phases, giving larger EDMs than those in the SM. Some mismatches 
with observations are # 13 and the neutrino mass square ratio. However, they may not be 
serious since we can improve the situation by extending to incorporate also type II seesaw, 
preserving the minimal SO(10) model. Leptogenesis is discussed in section 2.6. In section 2.7 
we go further to the analyses of the Higgs superpotential and reveal the concrete structure 
of intermediate energy scales between GUT and the SM, which leads us to the ruin of gauge 
coupling unigfication of the MSSSM. Another deficit of minimal SO(10), fast proton decay, is 
discussed in section 2.8. So we need some modifications of the minimal SO(10) model. Some 
solutions are discussed in section 2.9 but other bad news of GUT in 4D (No-Go theorem 
on SUSY breaking) is explained. Consequently, in the susequent part, we consider a class 
of SO(10) models with 5D orbifold [p3l . In Chapter 3 we consider an SO(10) model in 5D 
and will explain how this model rescues the problems of minimal SO(10) GUT in 4D. In 
section 3.2, we explain the setup of our model, [Q, where all matters and Higgs multiplets 
reside only on a Pati-Salam brane (PS brane) where the PS gauge symmetry is manifest, so 
that a low energy effective description of this model is nothing but the PS model in 4D with 
a special set of matter and Higgs multiplets. At energies higher than the compactification 
scale, the Kaluza-Klein (KK) modes of the bulk SO(10) gauge multiplet are involved in the 
particle content. In section 3.3, the gauge coupling unification is shown to be successfully 
realized by incorporating the KK mode threshold corrections into the gauge coupling running. 
The unification scale (Mqut) and the compactification scale (M c ) which was set to be the 
same as the PS symmetry breaking scale (vps) is found to be Mqut = 4.6 x 10 17 GeV and 
M c = vps = 1-2 x 10 16 GeV. The improvements of mass spectra data fitting will be discussed. 
This is rather trivial fact. In section 3.4, we apply this SO(10) model to the inflationary 
|35fl . The idea of inflation []36| has been strongly favored from the view point of 



scenario 



not only providing the solutions to the horizon and flatness problems of the standard big 
bang cosmology but also recent precise cosmological observations on the cosmic microwave 
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background radiation and the large scale structure in the Universe. Therefore, it is an 
important task to construct a realistic inflation model based on some well-motivated particle 
physics model. Single-field inflation is disfavored by its requirement of tiny parameter value 
to reproduce the results of COBE |37| and WMAP Hybrid inflation @ [0 |g| solves 
the above problem but gives rise to monopole problem in the original form. Variants of hybrid 
inflation, in particular, applicable to SUSY GUT models have been proposed: standard j42| . 
shifted [43] and smooth [45[ hybrid inflation models. Some of these models are based on the 
SUSY PS model with one singlet and Higgs multiplets (4, 1, 2) © (4, 1, 2) whose VEVs break 
the PS symmetry to the SM one. Interestingly, except for the singlet field, the orbifold GUT 
model of Ref. [Q, which we are interested in, has the same particle content. Therefore, the 
GUT model can naturally incorporate hybrid inflation and is constrained from the inflation 
model. So far we have assumed for simplicity M c = vps- However in this case LSP becomes 
stau. In section 3.5, we show that the neutralino becomes the LSP by generalizing M c ^ vps 
and DM candidate 
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In sections 3.6 and 3.7 we reanalyze leptogenesis |47j and LFV, 
which were discussed in the minimal SO(10) model in 4D, respectively in the scheme of 
SO(10) in 5D. According to the recent discovery of Higgs (like) particle, we have added a 
subsection on the impact of the LHC results to GUT. 

The last section is devoted to discussion. Appendix is served for a compact mathematical 
resume of SO(10) group property. 
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Chapter 2 

Renormalizable Minimal SO(10) GUT 
in Four Dimensions 



2.1 Why Do We Need GUT ? 

It is natural to start our review with this title since GUT is not necessarily indispensable 
for new physics BSM to all model builders. In going beyond the SM, we have rather serious 
constraints on matter contents, whereas we have no definite criteria for Higgs sector. One 



of the reasons for more Higgs than the SM is concerned with electro- weak baryogenesis [48 
We need strong first-order phase transition at the sphaleron transition (its rate T), 



It requires in the SM 



2E 



> 1 



(2-1) 



(2.2) 



with 



E 



2nv 3 



(6m 3 w + 3m|) 



(2.3) 



and A is the coefficient of ^0 4 . This leads us to m# < 56 GeV. If we add additional scalar 



particles, E can get larger and m# can go beyond the LEP or LHC bound [49 



So it is rather natural to consider more Higgs than one SM Higgs. In this case we have 
two major constraints: One is the p parameter, 



P 



m 



m z cofi 2 6w 



(2.4) 



and the other is FCNC. In two and more Higgs doublet model, FCNC appears at tree level 
in general. Let us consider p parameter first. It is famous that this parameter first predicted 
correctly large top quark mass before its discovery. Indeed, one loop quark correction to W 
propagator gives pUJ 



p = l + 



3G f 



8?r 2 \/2 



2 i 2 

m t + m b 



2m 2 t m 2 b trr^ 

2 2 1 2 

m t ~ m b \ m b 



2\ 1 



(2.5) 
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Two Higgs Doublet Model (2HDM) [|51[] contributes at one loop level and Higgs Triplet 



Model (HTM) @ f5| does at tree level. 

So first we consider on the HTM. The HTM is the minimum extension to the SM. It adds 
only one Higgs triplet and no matter field even right-handed neutrino. In the Higgs Triplet 
Model, we introduce a SU(2) triplet Y = 2 scalar as 

/ A + /\/2 A ++ \ 2 

A = ( _ A + ) ' ^ tri P let Yukawa = —hapL^ia APj^L^ + h.c. (2.6) 

This model generates neutrino masses without right-handed neutrinos with the triplet vac- 
uum expectation value v a which is given by the explicit breaking of the lepton number. This 
model is very predictive because of a clear relation 



V2v A h aP , (2.7) 



where m a p denotes the Majorana mass matrix for neutrinos. 
The experimental limit of p = 1.0004+£;^ at 2a f54[ gives 



— < 0.01, (2.8) 
vh 

where v = 246 GeV. 

Next simple model is two Higgs doublet model (2HDM), which add, in addition to the 
SM Higgs doublet Hi, another Higgs doublet H 2 . There are several types of the model 
depending on which doublet couples with which fermion: 

type I (SM-like) : Hi couples with all fermions 
H 2 decouples with fermions 
type II (MSSM-like) : Hi couples with down-type quarks and charged leptons 

H 2 couples with up-type quarks 
type III (general) : both of Higgs doublets couple with all fermions 
etc. etc. 



In the 2HDM as a whole, there is no criteria to classify matter content and it is the reason 
why there are several types mentioned above. Its strategy seems to let observations tell a 
story without model prejudice. In other words, in the 2HDM, it is rather difficult to specify 
the model and predict new phenomena. Also the masses of neutral and charged Higgses 

and phases are tightly constrained from Rb = T (z^hadrons) ' - > s l)i B — B mixing, p 
parameter etc., and we should take those constraints all into account. 

Both the HTM and 2HDM are very simple and useful towards the final theory. However, 
they are phenomenological models and can not be comprehensive new theories BSM. 
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2.2 What and Why Is SO(10) Group ? 



Chiral (left-handed) fermions in the SM with right-handed vr are composed of left-handed 

SU (2) l doublets (Q and L) and the charge conjugates of right-handed singlets ((ur) c , (cLr) , (e#) c , (vr) 

(hereafter we abbreviate them u c , .. for simplicty), 



Q 



u 



U r Uy Uf, 

y 

(u rc , u yc , u bc ) 



df dm d/jj 



3 ' 2 '« 

o 



3,1,— 
' ' 3 



d c = (d rc , d yc , d bc ) 



3,1,- 
' '3 



(2.9) 



(1,1,1), u c = (1,1,0). 



su(2) L x su(2) R m 



The first unfication scheme began with partial unification of quark- leptons under SU(A)c x 

(2.10) 



(4,2,1) 



Uj- Uy Ub V e 

d r d y db e 



Likewise, (4,1,2) is the charge conjugation of their right-handed partners. Here, L(R) 
indicates left-handed (right-handed) fermions and 1 indicates the first family. We have, of 
course, the second, third families. Thus lepton number was considered as fourth color. 
One of the great aqchievments of the PS model is the realization of charge quantization 

via 



Q = T L3 + ^- = T J 



L3 



T, 



B-L 



R3 



2 2 (2.11) 

SU(5) model unifies strong, electromagnetic, and weak forces but does not in single multiplet 
for matters matters, 

5 = (d rc , d yc , d bc } e- 



(2.12) 



and 
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1 

72 








u cr 


-u cy 


-u r 


-d r 


\ 




_ u cb 





u cr 


-u y 


-d y 






u cy 


-u cr 





-u b 


-d b 






U r 


u v 


u b 





-e+ 




V 


cT 


d» 


d b 


e+ 





/ 



(2.13) 



and 



Minimal SU(5) model gives problematic mass relation 

M e = M d 



(2.14) 
(2.15) 
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and fast proton decay. The flipped SU(5) JL4| circumvents this pathology and gives a good 
instrument, missing partner mechanism, for doublet-triplet problem (See subsection 2.9.2). 
Unfortunately, v c moves into 10-plet from SU(5) singlet and heavy Majorana mass term 
appears as unrenormalizable term. 

On the other hand, the fundamental representation 16 in SO (10) is an anomaly free and 
includes all fermions in a single multiplet 

16 L = (QX, cf , L, i/, e c ) L . (2.16) 

16 is decomposed as 

16 L = (4,2,1) + ((4,1,2) (2.17) 

under PS and as 

16 L = 5 + 10 + 1 (2.18) 

under SU(5), respectively. Of course, we need another Higgs to cancell the contribution of 
the fermion partners of the Higgs of SM if SUSY is involved. 

Let us go to further extension of gauge group. Next large gauge group is of rank 6, Eq, 
whose most conventional mass assignment is 



Here 

-U). *■-(£)' ^ 

and g is colored SU(2) singlet. In ( |2.19| ) the contents of upper line are matters of 16 = 
10 + 5 + 1 and those of lower are exotics of 10 = 5 + 5 and 1, respectively. N and iV c play 
the role of Higgs fields, H d and H u , respectively. 

Unlike the case of SO(10), it includes exotic particles. One of advantages was that Higgs 
fields are assigned in matter multiplet as H and H c (H implies the scalar partner of H). 
However, we must have three copies of Higgs multiplets corresponding to three generations. 
So some one consider that only and have vev and the other Hi (i = 1,2) and 
if? (i = 1,2) are non-Higgs ||57|| . So we must consider an additional scheme for the mass 
generation of the remaining first and second families. Also, if we set fermions and osons in 
the same multiplet we can not define R-parity, 

R = (_ 1 )3(B-D+25_ (2 21) 



R-parity is Z 2 subgroup of U(1)b-l which is carried out by 126 for SO(10) and 351 s = 
l(-8) + 10(-2) + 16(-5) + 54(4) + 126(-2) + 144(1) for E 6 . 

Also general E e invariant Yukawa coupling 27 x 27 x 27 leads us to B-L violation at low 
energy scale. To circumvent the troubles, one method is to prepare 27 Higgs additionally 
and separate matters and Higgs sectors. In that case we must add another Higgs to realize 
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non diagonal mixing angle of lepton mass matrix. For the renormalized case it may be 351//. 
As for unrenormalized effective action we do not adopt such strategy as we have repeated 
several times. Anyhow, this is far from being minimal mentioned above, and hereafter we 
concentrate on SO(10) GUT. 

First, we give a brief review of the minimal SUSY SO(10) model. 

There are many excellent textbooks for GUT. Here we list some of them emphasizing 
neutrino physics, [EHJ and 
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2.3 Yukawa Coupling 

A major part of great successes of the SM comes from its perturbative predictions. From the 
bottom-up approach, we have seen in the introduction that theoretical developments overlap 
with the extension of renormalizable area. So at least the main part of GUT should be 
guided from renormalizability. If we respect renormalizability, the dimension of the coupling 
constant (we symbolically write it Y here) must have zero or negative index in length units. 
So if fermion masses are generated by Higgs mechanism as in the SM, their interaction must 
be Yukawa coupling. 

Y-fuiHuQj + Y^d t H d Q 3 + 
+Y i - lj) V i H u L J + V,}'T, //,,/., + h.c. (2.22) 

We have written here for the case of MSSM . In the case of the SM+massive (Dirac) neutrino, 
we may change 

F u ->$ = (0°,</r) T , H d ^ir 2 <$>*. (2.23) 
So matter multiplets products in Yukawa coupling become group theoretically, 

16® 16 = 10© 120© 126. (2.24) 



The details of group theoretical arguments of SO(10) are given in |6(| |6T[ 



So the Higgs fields which can construct SO (10) singlet with bi-product of above fermions 
are 10, 120, and 126. Obviously single Higgs is incompatible with the observed mixing ma- 
trices, CKM and MNS. Since 10 is the simplest and inevitable, so the set of {10 and 120} or 
{10 and 126} is a minimal model. 126 is more essential than 120} in connection with neu- 
trino mass since the former includes (10, 3, 1) and (10, 1, 3) under SU(A) C ®SU (2) L ®SU{2) R . 
These two subgroups, if they have vevs, induce the right-handed Majorana and left-handed 
Majorana neutrino, respectively. It also automatically conserves R-parity f6"2"f . (Whereas 
SU(5) GUT also induces R-parity violation term via 105 5 term like LLe c , LQd c , u c d c d c .) 

So we select {10 and 126} in the Yukawa coupling. This model is called the renormal- 
izable minimal SUSY SO(10) GUT (the minimal SO(10) GUT). 

This model was first applied to neutrino oscillation in ||16|| . However, it did not reproduce 
the large mixing angles. It has been pointed out that CP-phases in the Yukawa sector 



play important roles to reproduce the neutrino oscillation data [17|. More detailed analysis 



incorporating the renormalization group (RG) effects in the context of MSSM has explicitly 
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shown that the model is consistent with the neutrino oscillation data at that time, Thus the 



minimal SO(10) GUT became a realistic model [18]. We give a brief review of the minimal 

SO(10) model. 

Yukawa coupling is given by 

W Y = yglGitfioie,- + y i i il6 i £Ti a6 16 i , (2.25) 

where 16j is the matter multiplet of the z-th generation, H w and Hi 2 q are the Higgs multiplet 
of 10 and 126 representations under SO(10), respectively. Note that, by virtue of the 
gauge symmetry, the Yukawa couplings, Y w and Y\i%, are, in general, complex symmetric 
3x3 matrices. After the symmetry breaking of SO(10) to SU(3) C x SU(2) L x U(l)y via 
SU(4) C x SU(2) L x SU(2) K or SU(5) x U(l), we find that two pairs of Higgs doublets in the 
same representation appear as the pair in the MSSM. One pair comes from (1,2,2) C 10 
and the other comes from (15, 2, 2) C 126. Using these two pairs of the Higgs doublets, the 
Yukawa couplings of Eq. (|2.25|) are rewritten as 



W Y = {U% (Y«H? + Y^H^) Q + (D% (Y^Hf + Y^H d 26 ) Q j 
+ (N^ (ygfl* - 3YH 6 H? 26 ) Lj + {E% {Y^H d 10 - SV&Ht*) Lj 
+ U (Yii 6 v L ) Lj + (N^ (y& v R ) (N c )j . (2.26) 

Here ur, cIr, u r and e# are the right-handed SU(2) L singlet quark and lepton superfields, Q 
and L are the left-handed SU(2)x doublet quark and lepton superfields, and H^t are 
up-type and down-type Higgs doublet superfields originated from H 10 and H U 6, respectively, 
and the last two terms are the Majorana mass terms of the left- and right-handed neutrinos 
developed by the VEVs (vl and vr) of the (10, 3, 1) and (10, 1, 3) Higgs. The factor —3 in 
the lepton sector is the Clebsch-Gordan coefficient. 

In order to preserve a successful gauge coupling unification, suppose that one pair of 
Higgs doublets given by a linear combination H^ d and is light while the other pair is 
heavy (> M GU t)- The light Higgs doublets are identified as the MSSM Higgs doublets (H u 
and H d ) and given by 

H d = a d H d + p d H d 2fi , (2.27) 

where a Utd and Utd denote elements of the unitary matrix which rotate the flavor basis in 
the original model into the (SUSY) mass eigenstates (see (|2.123|) in detail). Omitting the 



heavy Higgs mass eigenstates, the low energy superpotential is described by only the light 
Higgs doublets H u and H d such that 

W Y = (U^ {a u YH + /rYj 6 ) H u Qj + {D% (a d Y? + (3 d Y^) H d Q j 

+ (N% {a u Y^ - ?>(5 u YH b ) H u Lj + (E% (a d Y^ - 3/3 d Y^) H d L 3 (2.28) 
+ Li (yj 6 v L ) Lj + (N^ (Y&v R ) (N c )j , 

where the formulas of the inverse unitary transformation of Eq. (|2.27|) , H^ d = a u,d H U)d + • • • 
and H±2Q = f3 u,d H U)d + ■ • • , have been used. Note that the elements of the unitary matrix, 
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a u ' d and f3 u ' d , are in general complex parameters, through which CP-violating phases are 
introduced into the fermion mass matrices. 

Providing the Higgs VEVs, H u = vsm(3 and H d = v cos (3 with v = 174GeV, the quark 
and lepton mass matrices can be read off as 

M u = C W M W + Ci2 6 Mi 26 , 

M d = M 10 + M 126 , 

M D = c 10 M 10 - 3c 126 M 126 , (2.29) 
M e = M 10 -3M 126 , 
M L = c L M 126 , 
M R = c R M 126 . 

Here M u , M d , M D , M e , M L , and M R denote the mass matrices of up-type quark, down-type 
quark, Dirac neutrino, charged-lepton, left-handed Majorana, and right-handed Majorana 
neutrino, respectively. Note that all the quark and lepton mass matrices are characterized 
by only two basic mass matrices, M 10 and Mi 26 , and four complex coefficients cio, ci 2 6, cl 
and Cr, which are defined as M w = Y w a d v cos (3, Mi 26 = Yyi%fi d v cos 0, cio = (a u /a d ) tan/3, 
c i26 = (/3 u //3 d ) tan/3, cl = vl/ ((3 d v cos (3)) and c R = v R /((3 d v cos/3)), respectively. These are 
the mass matrix relations required by the minimal SO (10) model. 

You should remark especially the relation between M d and M e different from (|2.15 ). 
In other word, Yio and Y\i§ must be of the same order at least in the first 2x2 matrix 
components. Another essential point is that the heavy right-handed is represented by M\2§ 
(also a part of the mass matrices of charged fermions). This is the main reason of its high 
predictivity of this theory. 

In the following in Part I, we set cl = as the first approximation. Except for c R , which 
is used to determine the overall neutrino mass scale, this system has fourteen free parameters 
in total if we resrict ourselves Mw and Mi 2 6 are real JT7|. This few parameters in all mass 
matrices assures the strong predictability of the minimal SO(10) Model. 

2.4 Data Fitting 

Data fitting is performed as follows: Firstly we fit the data of charged fermions, masses of 
quarks and charged leptons and CKM mixing angles. Using the parameters fixed by this 
process we proceed to fit the neutrino data. 

2.4.1 Charged fermions 

Eliminating Mio and M\2q from Eq.( |2.29D , we obtain 



M e = c d M d + c u M u , (2.30) 

where 

3Cio + Ci 26 4 

Cd = , c u = . (2.31) 

ClO — C126 Cio — c 126 
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Since M u , M d , and M e are complex symmetric matrices, they are diagonalized by unitary 
matrices U u , U d , and U e , respectively, as 



U^M U U U = D U , UjM d U d = D d , UjM e U e = D ei 



(2.32) 



where D u , D d , and D e are diagonal matrices given by 



D u = diag(m u , m c , m t ) , D d = diag(m d , m s , m b ) 
D e = diag(m e , m^m T ) . 



Since the CKM matrix V q is given by 



V n 



u'Em 



the relation ( [2.30 ) is re- written as follows: 

(UlU u ) T D e (UlU u ) = c d V q D d V q T + c u D u . 
Therefore, we obtain the independent three equations: 

TrD e Dt = \c d \ 2 Tr[(V q D d V^ + KD u )(V q D d V^ + KD u y 
Tr(D e Dl) 2 = \c d \ 4 Tr[((V q D d V q T + nD u )(V q D d V q T + kD^) 2 
deW e D\ = M 6 det \(V q D d V q T + nD u )(V q D d V q T + kD u )^ 



(2.33) 
(2.34) 
(2.35) 

(2.36) 
(2.37) 
(2.38) 



with k = c u /c d . By eliminating the parameter c d , we have two equations for the parameter 



(ml + m 2 + m 2 .) 3 

(ml + m 2 + ml) 2 
2(mlm 2 + m 2 m 2 T + m 2 T m 2 ^) 



(12^61) 3 

(EH) 



(12^61) 2 - ( p7|) 



(2.39) 
(2.40) 



Here ( |2.36| ) 3 , for instance, means the right-hand side of Eq. Q2.36D to the third power. Let 
us denote the parameter values of k evaluated from Eqs.( p.39"|) and ( p.40|) as ka and Kb, 
respectively. If ka and Kb coincide with each other, then we have a possibility that the 
SO(10) GUT model can reproduce the observed quark and charged lepton mass spectra. If 
ka and Kb do not so, the SO(10) model with one 10 and one 126 Higgs scalars is ruled out, 
and we must bring more Higgs scalars into the model. The result is depicted in Fig.|2~l. 



2.4.2 The number of parameters in the minimal SO(10) model 

As we have discussed in the previous section, among four freedoms of complex {c w , Ci 2 e} or 
{c d , k}, we have been able to fix the three of them, k and \c d \. This is not accidental. Let 
us discuss the situation in details in the SO(10) two Higgs model. 
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Figure 2.1: Contour plot on complex K-plane. The vertical line and the circle correspond 
to the solutions of Eqs. (|2.39|) and ( |2.40|) , respectively. 



By using the relation ( [2.35| ), we have investigated whether there is a set of parameters 
which can give the 13 observable quantities D e , D u , D d , and V q or not. We can rewrite 
Eq.(p35|) as 

A T e D e A e = c d (V q D d V q T + kD u ), (2.41) 

where 

A e = Up u , (2.42) 

c d = \c d \e i(T . (2.43) 

The quantities D e , D u , D d , and V q are inputs, and the quantities \c d \, k, and A e are the 
parameters which should be fixed from those observed quantities. In general, an n x n unitary 
matrix for n generations has n 2 parameters. Therefore, the number of the parameters is 

7V(pmt) = N(A e ) + N(c d ) + N(k) = n 2 + 2 + 2. (2.44) 

On the other hand, the number of equations is 

iV(eqs) =n(n+l), (2.45) 

because Eq.( ^.41 ) is symmetric. Therefore, the number of the unfixed parameters is given 
by 

7V free = 7V(pmt) - iV(eqs) = 4 - n = 1, (2.46) 
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for n — 3, i.e., the 13 observed quantities fix the parameters \cd\, K, and A e , but 1 parameter 
a remains as an unknown parameter [F7] and [T$\. 

Let us reconsider this number counting in more detail, going back to (|2.29|) . We first set 
the base as M u = D u (3). M d is 3 x 3 complex symmetric matrix (12) which is diagonalized 
by a general unitary matrix V as 

M d = V*D d Vl (2.47) 
Here nine parameters of V are divided into the CKM matrix (4) and five phases, 

V = e ia e if3T3 e^ Ts V CKM e^' T3 e^' n . (2.48) 

However, the first three phases are delieted by rephasing. Thus we have 3+6=9 in M d . 
Furthermore we have 2 complex number c d , k (4). Thus totally 3+9+4=16. In this numbers 
we have set the above two phases (/?', 7') as or 7r for simplicity. Thus we have 14 parameters. 
[] Using the remaining one parameter and one cr (in ( |2.29 )) we can fit all mass matrices 



including Mu and Mr, that is, all low energy phenomena. This is a miraculous predictivity. 

2.4.3 Application to neutrino sector 

Next, we proceed to study how to predict neutrino masses 

D v = UjM u Uu (2.49) 

and the Maki-Nakagawa-Sakata (MNS) mixing matrix 

U e = UjU* u = A* e Al (2.50) 

by using the observed quantities D e , D u , D d , and V q and the parameter values \c d \, k, and 
A e fixed by Eq.(^4Tp. 



SO(10) GUT asserts that the Dirac neutrino mass matrix is given by the form 

M D = c 10 M 10 - 3c 126 M 126 (2.51) 

and Majorana mass matrices of the left-handed and right-handed neutrinos, Ml and Mr, 
are proportional to the matrix M\2§- 

M L = c L M 126 , M R = c R M 126 , (2.52) 

where M\q and M\2q are related to the quark and charged lepton mass matrices M u , M d , 
and M e as follows: 

Ml0 = 3 -M±±^, (2.53) 

Mue = Mi^k, (2 . 54) 



1 This is in contrast with 28 parameters in the SM+iA 
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Table 2.1: The input values of tan (3, m s (Mz) and 5 in the CKM matrix and the outputs for 



the neutrino oscillation parameters. 



tan/3 


m s (M z ) 


5 


a 


sin 2 29 12 


sin 2 2fl 23 


sin 2 2#i3 




40 


0.0718 


93.6° 


3.190 


0.738 


0.900 


0.163 


0.205 


45 


0.0729 


86.4° 


3.198 


0.723 


0.895 


0.164 


0.188 


50 


0.0747 


77.4° 


3.200 


0.683 


0.901 


0.164 


0.200 


55 


0.0800 


57.6° 


3.201 


0.638 


0.878 


0.152 


0.198 



Then the neutrino mass matrix derived form the seesaw mechanism becomes 

M u = M L -M D M^M T D 
= c L M 126 



-c 



R 



\c w M w - 3c 126 Mi 26 )M^(c 10 M 1 o - 3c 126 M 126 ) T . (2.55) 



In the present paper we adopt cl = 0. Also we may ignore the phase of cr which does not 
affect the observed values. Therefore, we can rewrite Eq. fl2.55|) as 



similarly to Eq. Q2.41 ), where 



with 



\c R \A T v D v A v = M D M^ l Ml, (2.56) 



M D = doM - 3ci26Mi, ( 2 - 57 ) 
M 10 = ^(3M d + M e ), (2.58) 

M 126 = i(M,-M e ), (2.59) 



M d = U^M d U u = V g D d V q T , (2.60) 



M e = UlM e U u = A T e D e A e 

= c d {V q D d V q T + kD u ). (2.61) 

The reasonable results we found at that time are listed in Table ^.1| . The shortening of almost 
all the minimal SO(10) models is large value of (9 13 (see Eq. (|1.3j ) for the recent result). 

As mentioned above, our resultant neutrino oscillation parameters are sensitive to all 
the input parameters. In other words, if we use the neutrino oscillation data as the input 
parameters, the other input, for example, the CP-phase in the CKM matrix can be regarded 
as the prediction of our model. It is a very interesting observation that the CP-phases listed 
above are in the region consistent with experiments. The CP-violation in the lepton sector 
is characterized by the Jarlskog parameter defined as 

J C P = lm[U e2 U; 2 U: 3 U, 3 ] , (2.62) 

where Ufi is the MNS matrix element. 
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Figure 2.2: Three mixing angles in the MNS matrix as functions of er[rad]. The graphs 
with the highest, middle and lowest peaks correspond to sin 2 2# 2 3, sin 2 2# 12 and sin 2 2# 13 , 
respectively. The plots of sin 2 2# 23 and sin 2 2# 13 have the sharp peaks at cr ~ 3.2[rad], while 
sin 2 2#i2 has the sharp peak at a ~ 3.3[rad] cited from JH3]. 



2.5 Lepton Flavour Violation and Dipole Moments 

Lepton flavour violation (LFV), anomalous magnetic dipole moment (anomalous MDM), 
and electric dipole moment (EDM) are discussed in the unified way. These phenomena are 
very sensitive to New Physics BSM via new CP phases and new particles. 

The SM gives negligibly small LFV probability in charged leptons, even taking into 
account the neutrino oscillation, 

Br(l -> oc ( ^] < 10" 54 , (2.63) 
\ m wJ 

and observation of LFV process becomes a clear signature of BSM. 

At tree level in the SM, the interaction of fermion ip (of mass and electro-magnetic 
charge eQ^) with photon is given by 

-eQ^YVA, = -^Hid^-i^^A^-^^a^F^ (2.64) 



2m^ 2m^ \ 3 ■ B 

It is clear in eq. (|2.65|) that a fermion has a MDM with g = 2 at tree level in the SM. 

In loop level in the SM and/or models BSM, following effective interaction of gauge 



19 



invariant form can be obtained: 



-i^ (A l lP L + A%P R ) a^ 3 F, u 



-j{A$ + A%)^a^F, u + - Ai)^^F, v 

= (A*i+Aw(* B ff °g)^ + i(A%-A^^ ff ° s ff Q E )^ (2-66) 

For the EDM and MDM, we take zero momentum of the photon. Then imaginary part of 
coefficients of the effective interaction vanishes because of the optical theorem (imaginary 
part of the forward scattering amplitude is given by the sum of possible cuts of intermediate 
states). We have an anomalous MDM and EDM d^ as 

H = 9 ^ = ~^MA% + A% (2.67) 



dj, = 21m(A%-Af). (2.68) 

Note that A^ and A R must include a fermion mass (m^ or fermion mass in the loop) because 
the effective interaction ipa^ip changes the chirality which can be done by the mass term 
in the fundamental Lagrangian. If one of particles in the loop is much heavier than others, 
A L and A R are suppressed by the mass. Thus, for large A L and/or A R , it is preferred that 
masses of particles in the loop are similar to each other. 



The explicit formulas of Al,r etc. used in our analysis are summarized in |63| |p4 
According to these papers, hereafter we renormalize Al and A R as 

£eff = -i^m l ~l j a flv F flv (A j ^P L + A»P R )U. (2.69) 

The decay rate of the LFV of charged lepton is given by 

m in) = ^< (l^l 2 + 14?) , (2-70) 

while the real diagonal components of Al,r contribute to the anomalous MDMs and EDMs 
of the charged-leptons such as 

«f SY = ^ = -mlBe[At + A«A (2.71) 
4/e = -m h Im(Al - A%) (2.72) 

Let us consider first an anomalous MDM of muon, a M = (g — 2)/2. The muon anomalous 
MDM has been measured very precisely p5| as 

a^ xp = 11659208.0(6.3) x 10~ 10 , (2.73) 
where the number in parentheses shows la uncertainty. On the other hand, the SM predicts 

aJ M [r] = 11659193.2(5.2) x KT 10 , (2.74) 
a J M [ e+e_ ] = 11659177.7(5.1) x 10~ 10 , (2.75) 
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Table 2.2: The input values of tan/3 and the outputs for the CP-violating observables 



tan (3 


(m v ) ee (eV) 


Jcp 


e 




40 


0.00122 


0.00110 


7.39 x 10" 


-5 


45 


0.00118 


-0.00429 


6.80 x 10" 


-5 


50 


0.00119 


-0.00631 


6.50 x 10" 


-5 


55 


0.00117 


-0.00612 


11.2 x 10" 


-5 



where the hadronic contributions to o^ m [t] and a^ M [e + e~] were calculated by using data 
of hadronic r decay and e + e~ annihilation to hadrons, respectively (see also |J7], [)8|, |69|, |70], 
[7lH). The deviations of the SM predictions from the experimental result are given by 

Aa M [r] = a^ xp - af 1 [r] = 14.8(8.2) x 10~ 10 , (2.76) 
Aa^[e + e-] = a^ xp - af*[e + e-} = 30.3(8.1) x lO" 10 . (2.77) 

These values of Aa M [r] and Aa M [e + e~] correspond to 1.8a and 3.7a deviations from the 
SM predictions, respectively. In order to clarify the parameter dependence of the decay 
amplitude, we give here an approximate formula of the LFV decay rate ||63|| . 



m ^ Iff) ~ ^rn% x 







2 














16vr 2 


MJ 





tan 2 /3, (2.78) 



where Ms is the average slepton mass at the electroweak scale, and ( Am? ) is the slepton 



mass estimated in Eq. ( p. 87 ). We can see that the neutrino Dirac Yukawa coupling matrix 



plays the crucial role in calculations of the LFV processes. We use the neutrino Dirac Yukawa 
coupling matrix of Eq. (|2.85|) in our numerical calculations. 

These quantities are evaluated by using the outputs presented in Table pTT| and the 
results are listed in Table pT3j p||9~lfl. 



Here some comments on the rate of the LFV processes and the muon g — 2 are added. 
The evidence of the neutrino flavor mixing implies that the lepton flavor of each generation 
is not individually conserved. Therefore the lepton flavor violating (LFV) processes in the 
charged lepton sector such as \i — > ej, t —> are allowed. In simply extended models so as 
to incorporate massive neutrinos into the SM, the rate of the LFV processes is accompanied 
by a highly suppression factor, the ratio of neutrino mass to the weak boson mass, because 
of the GIM mechanism, and is far out of the reach of the experimental detection. However, 
in supersymmetric models, the situation is quite different. In this case, soft SUSY breaking 
parameters can be new LFV sources, and the rate of the LFV processes are suppressed by 
only the scale of the soft SUSY breaking parameters, which is assumed to be the electroweak 
scale. Thus the huge enhancement occurs compared to the previous case. In fact, the 
LFV processes can be one of the most important processes as the low-energy SUSY search. 
However, this needs another assumption of universal boundary condition of supersymmetry- 
breaking parameters independent on the GUT framework. 
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Universal SUSY-breaking 

Csoft = -~(M 3 gg + M 2 WW + M 1 BB + c.c. 



(usl u QH u - da d QH d - ea e LH d + c.c.j (2.79) 



Q^rrigQ — L^m 2 L L — urn^u — dm d d — era^e 

,2 TJ* TT „,2 



- m 2 Hu H* u H u -m 2 Hd H* d H d -{bH u H d + c.c). 

So generically we have 19 parameters (3 gaugino masses + tan/3 + fi + vn^ + 10 sfermion 
masses + 3 trilinear terms), called phenomenological MSSM (pMSSM). Universal SUSY- 
breaking is a very strong assumpotion that it requires not only flavour-blindness but also 
universality over quarks and leptons at GUT scale, 

m Q = m l = m | = m | = m i = ?7 ^ 1 3» ( 2 - 80 ) 

m Hu = m Hd = m (2.81) 



M 3 _ M 2 _ Mi _ M 1/2 



£ 2 2 0i ^ 



(2.82) 



a u = A Y u , ad = 4, Y d a u = A Y e . (2.83) 

This MSSM+universal SUSY breaking is called constrained MSSM (CMSSM) J72|. If, in 
place of ( |2.81|) , we set m Hu , m Hd as free parameters, it is called non- universal Higgs masses 



(NUHM2) model |73| or called NUHM1 for m Hu = ±m Hd WM. From the SO(10) view point, 



there is some reason to set the universality between squarks and sleptons at GUT scale 
since all matters belong to a single 16 but there is no definite reason to extend it to Higgs 
masses. (We will argue on the problems on CMSSM or its alternatives at the last part of 
this review in connection with Higgs-like boson around 125 GeV discovered at the LHC). 
Hereafter, we will discuss in the CMSSM framework. We evaluate the rate of the LFV 
processes in the minimal SUSY SO(10) model, where the neutrino Dirac Yukawa couplings 



are the primary LFV sources. Although in Ref. |L8| various cases with given tan j3 = 40 — 55 
have been analyzed, we consider only the case tan (3 = 45 in the following. Our final result is 
almost insensitive to tan (3 values in the above range. The predictions of the minimal SUSY 
SO(10) model necessary for the LFV processes are as follows (T^|: with o = 3.198 fixed, the 
right-handed Majorana neutrino mass eigenvalues are found to be (in GeV) 

M Rl = 1.64 x 10 11 , M R2 = 2.50 x 10 12 and M Rz = 8.22 x 10 12 , (2.84) 

where c R is fixed so that Am^ = 2 x 10~ 3 eV 2 . In the basis where both of the charged- 
lepton and right-handed Majorana neutrino mass matrices are diagonal with real and positive 
eigenvalues, the neutrino Dirac Yukawa coupling matrix at the GUT scale is found to be 

-0.000135 - 0.00273i 0.00113 + 0.0136z 0.0339 + 0.0580i 
Y v = | 0.00759 + 0.0119* -0.0270 - 0.00419z -0.272 - 0.175i J . (2.85) 
-0.0280 + 0.00397z 0.0635 - 0.0119? 0.491 - 0.526i 



2 



We are now reconsid ering data fitting with the update experimental data and new RGE results. It gives 



the differen values from fl2.85|) but the LFV results are not essentially changed. 
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LFV effect most directly emerges in the left-handed slepton mass matrix through the RGEs 
such as f6"3" | 



MSSM 

+ ^ {m}YX + Y}Y v m\ + 2YjmlY l/ + 



(2.86) 



where the first term in the right-hand side denotes the normal MSSM term with no LFV. We 
have found Y v explicitly and we can calculate LFV and related phenomena unambiguously 



18|j . In the leading-logarithmic approximation, the off-diagonal components (i ^ j) of the 



left-handed slepton mass matrix are estimated as 

where the distinct thresholds of the right-handed Majorana neutrinos are taken into account 
by the matrix L = log [M GUT 

The recent Wilkinson Microwave Anisotropy Probe (WMAP) satellite data [38| provide 



estimations of various cosmological parameters with greater accuracy. The current density 
of the universe is composed of about 73% of dark energy and 27% of matter. Most of the 
matter density is in the form of the CDM, and its density is estimated to be (in 2a range) 

n CD Mh 2 = 0.1126ljj;ffi • (2.88) 

The parameter space of the CMSSM which allows the neutralino relic density suitable for 
the cold dark matter has been recently re-examined in the light of the WMAP data [75|. It 



has been shown that the resultant parameter space is dramatically reduced into the narrow 
stripe due to the great accuracy of the WMAP data. It is interesting to combine this result 
with our analysis of the LFV processes and the muon g — 2. In the case relevant to our 
analysis, tan /3 = 45, fi > and A = 0, we can read off the approximate relation between 
m and M 1( / 2 such as (see Figure 1 in the second paper of Ref. 



m (GeV) = — M V2 (GeV) + 150(GeV) , (2.89) 
28 

along which the neutralino CDM is realized. M\/2 parameter space is constrained within the 
range 300GeV < Mi/ 2 < lOOOGeV due to the experimental bound on the SUSY contribution 
to the b — > sj branching ratio and the unwanted stau LSP parameter region. We show 
Br(/i — > ej) and the muon g — 2 as functions of Mi/ 2 in Fig. |2.3| along the neutralino CDM 
condition of Eq. ( |2.89p . We find the parameter region, 560GeV < Mi/ 2 < 800GeV, being 
consistent with all the experimental data. 

The semileptonic flavor violation processes were also considered in [pip , for instance, 



for r — > e (fi ) 7T°, r — > e (// ) r], r —> e (/j, )rf, t — > e (// ) p°, r — > e (fi ) (j), 
t~ — > e~(ji~) oj, etc. 
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Figure 2.3: The branching ratio 
2 in units of KT 10 



muon g 



5af USY 



Log 10 [Br(n — y ej)], b: the SUSY contribution to the 
\ , and c: the electron EDM, Log 10 [|d e |[e cm]]. 
All these figures are plotted as a function of M1/2 (GeV) along the cosmo logical constraint 
of Eq. ( |2.89 ). Trilinear term Aq is assumed to be zero except for the last panel. The last 
panel is added for the reference to see the behaviour of non zero Aq, where the branching 
ratios, Br(r — > ^7) (top) and Br(jU — > ej) (bottom) are given as functions of A (GeV) for 
m = 600 GeV and M 1/2 = 800 GeV. All are cited from PJ. 



Similarly we can estimate EDM from ( |2.68[ ) . Electron EDM is depicted also in Fig.[Q 
When the KamLAND data |23|] was released, the results in [|18l were found to be deviated 



by 3<t from the observations. Afterward this minimal SO(10) was modified by many authors, 
using the so-called type-II seesaw mechanism [76| and/or considering a 120 Higgs coupling 
to the matter in addition to the 126 Higgs [ff7| . Based on an elaborate input data scan [78 



79j, it has been shown that the minimal SO(10) is essentially consistent with low energy 
data of fermion masses and mixing angles. The importance of the threshold corrections 
was also discussed in [BU 



2.6 Leptogenesis 

Cosmic baryon asymmetry is one of the most important subjects for new physics BSM. 
Sakharov pointed out |y| that we need three conditions for baryogenesis: 



3 See ( |3.90| ) for the recent result of 
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1. Bayon number violation, 

2. CP violation, 

3. Out of equilibrium condition. 

As we emphasized in the previous section, CP violation process in the SM is parametrized 
by the Jarlskog parameter. Its magnitude is too small to generate the observed baryon 
asymmetry Y^b = 0(1CT 10 ). 

One of the reasons that we adopted 126 was that it includes (10, 1,3) which generates 
heavy right-handed Majorana mass term and induces A(B — L) = 2 with additional CP- 
violating (CPV) phases. 

The minimal SO (10) GUT has many scalar fields and many CPV phases, heavy right- 
handed neutrino Nr is not the unique parent for baryon asymmetry []. However, it seems to 
be most natural to accept leptogenesis via Nr |$3| . Supersymmetry requires the leptogenesis 



by the decays of both the lightest heavy right-handed neutrino Y^ f and sneutrino Yl b equal 
and total lepton asymmetry is f& 



Y L = Y Lf +Y Ls . (2.90) 
The processes of neutrino and sneutrino decays are essentially same and sneutrino case will 



be discussed at section 3.6. 



Thermal leptogenesis Nr has a large mass and no gauge interaction, and is out of 
equilibrium. The effective Lagrangian at energies lower than the right-handed neutrino 
masses is 



eff 



- J d 2 9 [y^N^H u + i M Ri NtNt\ + h. 



c. 



(2.91) 



where i,j = 1,2,3 denote the generation indices and Y v the Yukawa coupling, L and H u 
are the lepton and the Higgs doublets chiral supermultiplets, respectively, and Mm. is the 
lepton- number- violating mass term of the right-handed neutrino N. The peculiar properties 
of the minimal SO(10) are that we can fix Y^ and Mm unambiguously from the low-energy 
phenomenologies of quarks and leptons. The lepton asymmetry in the universe is generated 
by CP- violating out-of-equilibrium decay of the heavy neutrinos, N — >■ £lH* and N — >■ £lH u . 
The leading contribution is given by the interference between the tree level and one-loop level 
decay amplitudes, and the CP-violating parameter is found to be |[5[ 

e = - * t , £ Im [<y v Yl%] {f(M 2 Rj /M R1 ) + 2g(M 2 Rj /M 2 R1 )} . 
o-K{Y v Yv)n j=2 3 

Here f(x) and g(x) correspond to the vertex and the wave function corrections, 



(2.92) 



9{x) 



1 + x)ln 



x 



X 



X 



2(1 -x) 



(2.93) 



indeed, there are alternative approaches 
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respectively, and both are reduced to ~ — ^ ^ ^' ^ ^^is approximation, £ 
becomes 

e = — t — V Im [(y/t)21 ^ . (2.94) 



Using the mass of neutrino via Type-I see-saw mechanism, M v = —YjM^Y^^u) 2 , e is 



further written as 86 



3 M m Im [(KM;yJ) n ] 
Wit(H u ) 2 (Y v Yj) n 
3 m u3 M R1 6 eS 



(2.95) 



In the minimal SO (10) model we have the definite form of Y v and estimate these values 
unambiguously. We have assumed that the lightest Ni decay dominantly contributes to the 
resultant lepton asymmetry. In fact, this is confirmed by numerical analysis in the case of 

Using the above e, generated Yb is described 



hierarchical right-handed neutrino masses J84 
as 



Y 



B 



9* 



-d 



(2.96) 



where g* ~ 100 is the effective degrees of freedom in the universe at T ~ Mm, and d < 1 is 
so-called the dilution factor. This factor parameterizes how the naively expected value Yb ~ 
e/ g* is reduced by washing-out processes. We can classify the washing-out processes into two 
cases with and without the external leg of the heavy right-handed neutrinos, respectively. 
The former includes the inverse-decay process and the lepton-number-violating scatterings 
mediated by the Higgs boson such as N+£ L <-> q~R+qL, where qi and are quark doublet 
and singlet, respectively. The latter case is the one induced by the effective dimension five 
interaction, 



C 



N 



- (YjM^Y v ).. {UE u fC-\L 3 H u ) 



(2.97) 



after integrating out the heavy right-handed neutrinos. This term is nothing but the one 
providing the see-saw mechanism . The importance of this interaction was discussed in |37j . 



where the interaction was shown to be necessary to avoid the false generation of the lepton 
asymmetry in thermal equilibrium. While numerical calculations |84| are necessary in 
order to evaluate the dilution factor precisely, Yb ~ c/g* roughly gives a correct answer, and 
the washing-out process is mostly not so effective. Note that this is the consequence from 



the current neutrino oscillation data as explained in [89 



These quantities are evaluated by using the outputs presented in Table P7T] , and the 

930- 

Unfortunately, 

the CP parameter e is too large to be consistent with the observed baryon asymmetry. In 
order to circumvent this trouble we made use of another pair of SU(2) doublets appearing in 



results are listed in Table \L3 

{m u ) ee is the averaged neutrino mass appearing in the Qv/3/3 process. 



26 



Table 2.3: The input values of tan/3 and the outputs for the CPV observables 



tan (3 


(m u ) ee (eV) 


Jcp 


e 




40 


0.00122 


0.00110 


7.39 x 10" 


-5 


45 


0.00118 


-0.00429 


6.80 x 10" 


-5 


50 


0.00119 


-0.00631 


6.50 x 10" 


-5 


55 


0.00117 


-0.00612 


11.2 x 10" 


-5 



the minimal SO(10) model. We solved the Boltzman equation and obtained the consistent 
Yb [p2| . However, in this case we need the extra Higgs other than those in the MSSM, 



which may raise the other problems. So it is deserved to consider an alternative solution 
to this overproduction. On the other hand the gravitino problem forces us low reheating 
temperature less than the mass of Mri- If we believe it, the above problem becomes fake 
since thermal Nr are not generated in the reheating era. So the minimal SO (10) model 
itself drives us to the other approaches such as non-thermal leptogenesis scenario |93| or the 
Affleck-Dine mechanism [|9"4j| . In the next section, we discuss on the non-thermal leptogenesis 
scenario in the minimal SO (10) model. 

Non-Thermal leptogenesis Now we turn to the discussions of the non-thermal leptoge- 
nesis scenario ||93|| . In the non-thermal leptogenesis scenario, the right-handed neutrinos are 
produced through the direct non-thermal decay of the inflaton. 

Here we give a concrete model to specify the inflaton. We add a singlet chiral supermul- 
tiplet which plays a role of inflaton S The interaction Lagrangian relevant to the inflaton 
and the right-handed neutrinos is given by 

C s = ~\j d 2 6 (^MjS 2 + KSN?N?j . (2.98) 

When inflaton gets a VEV, it gives rise to the Majorana masses for the right-handed neutrinos 



in addition to the VEV of (10, 1,3) in 126 under SU(A) PS x SU(2) L x SU(2) R 0, P 



However, the VEV (S) is posted around the GUT scale and A; is found to be 10 later, and 
this contribution gives a tiny correction to Mr. Also the first term in Eq. ( p. 98 ) dominates 
over the second, and is reduced to the chaotic inflationary model [|95] . 
In such a superpotential, the inflaton decay rate is given by 



I A, 



12 



T(S -> NiNt) ~ . (2.99) 

Then the consequently produced reheating temperature is obtained by 

/ 4S \ 1/4 

If the inflaton dominantly couples to N, the branching ratio of this decay process is, of 
course, BR ~ 1. Then the produced baryon asymmetry of the universe can be calculated by 
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using the following formula, 



s J 



-0.35 x 



n Nl 



x 



n Nl 



-0.35 x - BR(S -> N^i, 



Tr 
Mi 



x e 



(2.101) 



With the hierarchical mass spectra for the right-handed neutrinos, it can be approximated 

as 



s J 



-1.95 x 10 -10 x BR x 



T 



R 



M. 



m 



10 6 GeV J V Mi J V 0.065 eV 



m„ 3 



x S, 



eff 



(2.102) 



where <5 e g = Im [(1^ / M*l^ , )ii] / [^3(1^1^)11] denotes the effective value of the CP violating 
phase parameter relevant to the leptogenesis and it can be estimated as 5 e s = —0.166 in our 
model. As it can easily be seen that it is possible to produce the baryon asymmetry of the 
universe by using the reheating temperature as low as, Tr < 10 6 GeV. Hence, a very wide 
range of the gravitino mass can be allowed, m 3 / 2 > 10 6 MeV. The result of the detailed 
numerical calculation based on Eq. ([2.101) is shown in Figs. 2. 4 and 2.5 p6|. 
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Figure 2.4: The predicted baryon asymmetry of the universe Yb = ub/s as a function of 
the inflaton mass Mj [GeV] with the reheating temperature Tr = 10 6 GeV and BR = 1. The 
vertical line represents the kinematical cut for the inflaton to have enough energy to decay, 
E > 2Mrx, i.e., the right side of this line is allowed from the kinematics. Two horizontal lines 
represent the upper and the lower bounds on the observed value of the baryon asymmetry 
at 95 % C.L. The lepton asymmetry parameter e has been taken from fl8| . 
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As shown in Fig.2.4 that the predicted inflaton mass is heavier than the lightest right- 
handed neutrino (Mri = 1.6 x 10 11 GeV) in our model Hence the non-thermal leptogen- 



esis is well workable. But in model p7\ , you can see from Fig.2.5 that the calculated inflaton 
mass is lighter than the lightest right-handed neutrino mass (Mri = 2.7 x 10 13 GeV), and 
the non-thermal leptogenesis scenario is prohibited by the kinematics. We hasten to add 
that this conclusion is valid under the non-thermal leptogenesis under the gravity mediated 
SUSY breaking scenario. 

It can be read from Fig.2.4 that the observed value of the baryon asymmetry leads to 
the inflaton mass around Mi ~ 5 x 10 11 GeV. This corresponds to the coupling constant 
of the inflaton to the right-handed neutrinos as Aj ~ 10~ 8 . Such a small coupling indicates 
that the model can naturally fit into the chaotic inflationary model |J5| based on a minimal 
supersymmetric SO(10) model. 



2.7 Higgs Superpotential- 

Symmetry Breaking Flows from GUT to the SM 

On the other hand, it has been long expected to construct a concrete Higgs sector of the 
minimal SO(10) model. 

The simplest Higgs superpotential The simplest Higgs superpotential at the renor- 
malizable level is given by fl99f , ||100|| 

W = mi$ 2 + m 2 AA + m 3 H 2 + Ai$ 3 + A 2 $AA + \ 3 $AH + X^AH , (2.103) 
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where $ = 210, A = 126, A = 126 and H = 10. The interactions of 210, 126, 126 and 
10 lead to some complexities in decomposing the GUT representations to the MSSM and 
in getting the low energy mass spectra. Particularly, the CG coefficients corresponding to 
the decompositions of SO(10) — > SU(3)c x SU(2) L x U(l)y have to be found. This problem 
was first attacked by X. G. He and S. Meljanac ||101|| and further by D. G. Lee [|99] and by 



J. Sato ||102|| . But they did not present the explicit form of mass matrices for a variety of 
Higgs fields and also did not perform a formulation of the proton life time analysis. This 
is very laborious work and it is indispensable for the data fit of low energy physics. We 



completed this program in [£5j (see also f2q|, ||27| ). This construction gives some constraints 
among the vacuum expectation values (VEVs) of several Higgs multiplets, which gives rise 
to a trouble in the gauge coupling unification p§| . The trouble comes from the fact that the 
observed neutrino oscillation data suggests the right-handed neutrino mass around 10 13 ~ 14 
GeV, which is far below the GUT scale. Indeed ( |2.103| ) contains five directions which are 



singlets under SU(3)c x SU(2) L x U(l)y- Three of them are included in 210, 

k = (1234), 

2 = (5678 + 5690 + 7890), (2.104) 

03 = (12 + 34)(56 + 78 + 90). 



one in 126 and 



v R = (13579), (2.105) 



and one in 126 



^=(24680). (2.106) 
Due to the D-flatness condition the VEVs v r and Ur are equal, 

v R = vE. (2.107) 

This intermediate scale is provided by Higgs field VEV, and several Higgs multiplets are 
expected to have their masses around the intermediate scale and contribute to the running 
of the gauge couplings. 

We write down the VEV conditions which preserve supersymmetry, with respect to the 
directions 0i, 2 , 3 , and Vr, respectively. 

2^ + SA^+A^ = 0, (2.108) 



2 ^ 2 + 3Al l r t#J +A2 T^! = °' (2 ' 109) 

2^ + 3A, [ ^ + + X < = o, (2.110) 




™ 2 + A 2 (-^ + ^ + ^) = 0. (2.111) 

\ioVg 10V2 10 ' 
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Eliminating v R , 0i and 03 from Eqs. (T2.108D - fl2.lll ), one obtains a fourth-order equation 



in 02- The corresponding fourth-order polynomial in 2 factorizes into a linear and a cubic 
term in 2 . Linear term gives the solution of the fourth-order equation which is very simple, 
02 = —3-^/2^2/^2, but it preserves the SU(5) symmetry. Therefore, it is physically not 
interesting. The cubic term solutions lead to the true SU(3)c x SU{2) L x Z7(l)y symmetry. 
Here we consider only the solutions with \vr\ 7^ 0. Eliminating vr -vr, 0i and 02 from Eqs. 



fl2.108|) -( j2.111|) , one obtains a fourth-order equation in 



Us + ^) {8*1 - 15 M^l + UM\<t> 3 - 3M\ + (0 3 - M l ) 2 M 2 } = 0, (2.112) 

Mi = 12 ("^) > -Mi = 60 (^) ■ (2.H3) 



where 



Any solution of the cubic equation in 03 is accompanied by the solutions 

3 (Ml -5 <%) 
91 " v/6 (Mi - 0s) 2 ' 

1 (.Mi - 2^103 - 0|) 



x/2 (M- 



>3 



(2.114) 
(2.115) 



_ 5 f\A 03(^- 3 03)^ + 01) 

= SUJ (AWsF • (2 - 116) 



The linear term gives the solution of the fourth-order equation ( p,112| ) which is very sim- 



pie, 03 = -6 (^). It leads to 0! = -^6 (^), 2 = -3y/2 (^) and v^-^) 



(^) v ^ (w*a) _ ^ (^2)' som ^ on preserves the SU(5) symmetry. Therefore, it 

is physically not interesting. Then we proceed to the most important part of the SO(10) 
GUT. We can not show the details of the scenario but only show the essential part of it p5|. 



Would-be Nambu-Goldstone bosons At first, we list the quantum numbers of the 
would-be Nambu-Goldstone (NG) modes under SU(3) C x SU(2) L x U(1) Y . 

• [(3,2,|)©(3,2,-|)], 

• [(3,2,-|)©(3,2,|)], 
. [(3, 1,-1) ©(3, !,§)], 

• [(1,1, 1)0(1,1,-1)], 

• [(1,1,0)]. 

Total number of the NG degrees of freedom is: 12 + 12 + 6 + 2 + 1 = 33. The cubic term 
solutions lead to the true SU(3)c x SU{2) L x U(l)y symmetry. ( [2.1 16| ) gives heavy right- 
handed neutrino, and the coefficient of ( [2.29|) is also written in terms of 03. 
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Electroweak Higgs doublet In the standard picture of the electroweak symmetry break- 
ing, we have the Higgs doublets which give masses to the matter. These masses should be 
less than or equal to the electroweak scale. Since we approximate the electroweak scale as 
zero, we must impose a constraint that the mass matrix should have one zero eigenvalue. 
We define 



and 



Hl° - <S A„ EE A^g, A u EE Agjgj, $ u EE (2.117) 

Hf = < 2 %7 \ A, ee A[; 5 2 -|, A d ee AgJ"^, $ d ee $£ o 2 -|. (2.118) 



In the basis {H^°, A u , A n , the mass matrix is written as 



/ 2T7i3 -^302 A3 03 A402 A 4 </>3 XjVR 



M, 



doublet 



M<t>2 A403 



VW-lVK m 2 + ^%-^ 

A3<fe _ A3 03 °V A 2 02 1 A 2 03 A 2 ffl 



y to -to 2m 1 + to + to 

(2.119) 

The corresponding mass terms of the superpotential read 

W m = (Hl°,A u ,A u ,<S> u ) M doublet (H™,A d ,A d ,<S> d ) T . (2.120) 

The requirement of the existence of a zero mode leads to the following condition. 

det M doub iet = 0. (2.121) 

For instance, in case of A 3 = 0, m 2 + j|% — = 0, we obtain a special solution to Eq. 
( |2.121 ), while it keeps a desirable vacuum and it does not produce any additional massless 
fields. However, we proceed our arguments hereafter without using this special solution. 
We can diagonalize the mass matrix, Mdoubiet by a bi-unitary transformation. 

U* M doublet yt = diag ( , M 126 , M 2 , M 3 ). (2.122) 
Then the mass eigenstates are written as 



(H u , hi hi hi) = (^°,A U ,A U ,$ U ) U T , 

(H d , hi h 2 d , h*) = (tf d 10 ,A d ,A,,$ d ) V T . (2.123) 

Here H u , H d are MSSM light Higgs doublets. We get the explicit form of U and V from 
( |2.119 ), and thus we can connect the oscillation data with the GUT Yukawa coupling. Thus 



the intermediate energy scales are severely constrained from the low energy neutrino data, 
and the gauge coupling unification at the GUT scale may be spoiled. 

This fact has been explicitly shown in [p8 |, where the gauge couplings are not unified any 



more and even the SU(2) gauge coupling blows up below the GUT scale (Fig. 3). Thus the 



32 



Msusy 



Hps 

M PG M3211 M3221 M G 




detailed analyses of superpotential was the great progress but it reveals unambiguously the 
details of structure, which also uncovers pathologies. 

However, this is easily remedied by the addition of 120 Higgs in Yukawa coupling [|97|- We 
mean that the dominant part may be governed by the minimal SO (10) but such generalization 
does not spoil the renormalizable SO(10) GUT yet. 



2.8 Proton Decay 

One of the problems we encountered in the minimal SO(10) GUT is the fast proton decay. 
After the symmetry breaking from SO(10) to SU(3)c X SU (2) L x U(l)y, the generic Yukawa 
interactions between the matter fields and the color triplet Higgs fields are given by 

W Y = YHH T {q l i ] +u c i d c j )+YH^{q i i J +u c l d c j ) 
+ Yto H T Qftfc + u\e) + d\v\ 

+ Yii 6 A T (^q i q j + u<re c j + dy :j 
+ ^ 6 A;«e| + ^). (2.124) 

Here we have defined 

_ (3,1,1) (3,1,-i) _ =x (3,l,i) _ (3,1,-1) ^ = (3,1,-1) 

n T — n (6,1,1) ' —-"(6,1,1) ' T (6,1,1) ' — ^(6,1,1) ' T (10,1,3) ' {*-l*0) 

For later use we define 

A T = A (6,l,l) ' At = A (6,l,l) > A T = A (T0,l!3) > $ T = $ (15,l!3) > = $(15,1,3) ■ ( 2 - 126 ) 
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In the basis {Hjt, A^t, Ajt, A^}, the mass matrix reads 



/ 



^triple 



2m 3 



V 



V10 V30 

M<jn _i_ A 4 </> 2 



A3</>1 I A3 02 
) + 





i + 

10 V30 



f- 

lOv^ 



5^ 

where 7^44 = 2™! + ^ + ^ + ^ and m 55 = m 2 + ^ + |g§. 



A401 A402 

m 2 


_ Mvr_ 
lOv^ 

A2 03 
15^2 



v / 2A 4 



m 2 





A2</>3 



A2fj?. 

5^ 
^55 



/ 

(2.127) 



The corresponding mass terms of the superpotential read 

W m = (H T , A T , A T , %, A^) M trip , et (# T) A T , A T , $ T , A^ T . (2.128) 

Integrating out At, $t, and A^, we obtain the effective Yukawa interactions between the 
matter fields and the color triplet Higgs fields as 

W Y = Y^H T (q l £ J +u^)+Y^ e A T (q l £ J +u^) 



+ 



MO 



m 3 i 
m 33 
1 



2 

m 33 



Y& A T {u\e) + d<tf) 



(2.129) 



Then the effective mass terms for the remaining color triplet Higgs fields are written as 

W£ = H T (aH T + bA T ) 
+ A T (cH T + dA T ) 



(H T , A T ) M T ( g ) . 



(2.130) 



Eqs. ( p.l29|) and Q2.130| ) leads us to the effective dimension-five interactions after integrating 



out the remaining color triplet Higgs fields [|103 



W R 



lijkl , 



(2.131) 



inducing the dangerous proton decay. Here, Cl and Cr are given by the Yukawa coupling 
matrices at the GUT scale, Mq, 



•ykl 
1 126 



Cf\M G ) 



r 10 



m 13 
m 33 



1 126' 



m 32 
m 33 



1 126 



M^ 1 



Vkl 
J 126 



(2.132) 
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Note that 



Y 10 

^126 



Yu 
Yd 



Y u 
Y d 



(2.133) 



Thus we have 



{Yu> 1?) (AM T A T y 



ykl 

u 
Vkl 
1 d 



(2.134) 



The Yukawa coupling matrices, Y w and Yi 2 q, are related to the corresponding mass matrices 
M 10 and M 12 6 such that 



Yi 

Y\2% 



ClO 



a u v sin (3 

Cl26 

f3 u v sin (3 



-M 



10) 



M 



126, 



(2.135) 



with f ~ 174.1 [GeV]. Here a u and /3" are the Higgs doublet mixing parameters introduced 



in (|2.28| ), which are restricted in the range \a u \ 2 + |/3 U | 2 < 1. Although these parameters are 



irrelevant to fit the low energy experimental data of the fermion mass matrices, there are 
theoretical lower bound on them in order for the resultant Yukawa coupling constant not 
to exceed the perturbative regime. Since C10, C126, M w and Mi 2 6 are the functions of only 
a, we can completely determine the Yukawa coupling matrices once a, a u and (3 U are fixed. 
In order to obtain the most conservative values of the proton decay rate, we make a choice 
of the Yukawa coupling matrices as small as possible. In the following analysis, we restrict 
the region of the parameters in the range (a u ) 2 + ((3 U ) 2 = 1 (we assume a u and (3 U real for 
simplicity). Here we present examples of the Yukawa coupling matrices with fixed o = 71. 
For tan = 2.5 with a u = 0.031, we find 



Y 



10 



126 



-6 , 



0.000839 + 2.79 x 10 
0.00151 - 0.0000265 i 
0.000692- 0.000818 i 



-0.0000613 - 2.17 x HT r z 
-0.000111 + 1.94 x I0~ 6 i 
-0.0000508 + 0.0000600 i 



0.00151 - 0.0000265 i 0.000692 - 0.000818 i 
0.00479 + 0.0000811 i -0.0128 + 3.17 x 10~ 6 i 
-0.0128 + 3.17 x 10" 6 i 0.525 - 0.0420 i 

(2.136) 

-0.000111 + 1.94 x 10" 6 i -0.0000508 + 0.0000600 i 



and for tan/3 = 10 with a u = 0.111, 



-0.000428 - 2.46 x 10" 
0.000941 - 2.33 x 10 



0.000941 - 
1.42 x 10" 



0.00132 i 
(2.137) 



Y 



10 



0.00101 + 1.87 x I0~ b i 
0.00179 - 0.0000439 % 



0.00179 - 0.0000439 i 
0.00541 + 0.000141 % 



0.000348 - 0.00125 i -0.0154 + 5.21 x 10 



-6 , 



0.000348- 0.00125 i 
-0.0154 + 5.21 x 10~ 6 
0.530 - 0.0567i 



(2.138) 



35 



-0.000244- 5.21 x l(T 7 i -0.000436 + 0.0000107 i -0.0000847 + 0.000305 i 
r [2( , = [ -0.000436 + 0.0000107i -0.00164 - 0.0000154 i 0.00375 - 1.27 x 10 6 i 
-0.0000847 + 0.000305 i 0.00375 - 1.27 x 10~ 6 % -0.000684 + 0.00626 i 

(2.139) 

For the effective color triplet Higgsino mass matrix, we assume the eigenvalues being the 
GUT scale, Mq = 2 x 10 16 [GeV], which is necessary to keep the successful gauge coupling 
unification. Then, in general, we can parameterize the 2x2 mass matrix as 

M T = M G I 2 x U, (2.140) 

with the unitary matrix, 

U = e i V a a ( C0S9 Sm9 \ e V°*. (2.141) 

y — sm U cos v J ' 

Here we omit an over all phase since it is irrelevant to calculations of the proton decay rate. 
Now there are five free parameters in total involved in the coefficient C^ kl , namely, a, a u , 
8, ip and <p' . Once these parameters are fixed, C % [ kl is completely determined. 

The proton decay mode via the dimension five operator in Eq. (|2.131 ) with the Wino 



dressing diagram is found to be dominant, and leads to the proton decay process, p — > K + v. 
The decay rate for this process is approximately estimated as (in the leading order of the 
Cabibbo angle A ~ 0.22) 



X 



£(2311 _ ^1312 _|_ ^ /£2312 _ £<1322\ 
L L y L L ) 

x 5.0 x 10 31 [years" 1 /GeV]. (2.142) 
Here the first term denotes the phase factor and the hadronic factor, 

(3 H u L (k) = (0\d L u L u R \p(k)), (2.143) 
and H = 0.0096 [GeV 3 ] is given by lattice calculations p4| . A L ~ 0.32, A s ~ 0.93 are 



the long-distance and the short-distance renormalization factors about the coefficient C l [ k \ 
respectively. f\ The third term in the first line comes from the Wino dressing diagram, 
and ms is a typical sparticle mass scale multiplied by the ratio of a sfermion and Wino. 
In the case with the mass hierarchy between the sfermions and the Wino (fhf 3> M 2 ), we 
find ms ~ fhf x (fhf/M 2 ). In the following numerical analysis, we take fhf = 1 [TeV] and 
M 2 = 100 [GeV]. 

Now we perform numerical analysis. Note that because of the very constrained flavor 
structure of the minimal SUSY 5O(10) model we can give definite predictions for the proton 



As suggested in Ref. [ 105 1 , it might be proper to use the renormalization factors Al, As in [ 105 , 
which directly treats the renormalization of the Wilson coefficients itself. But here, we adopt the use of the 
conventional factors Al, As to compare our results to the previous ones. 
6 See the comments on the recent LHC results in the last section. 
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decay rate once the five parameters in the above are fixed. For a specific choice of the Yukawa 
coupling matrices in the minimal SO(10) model with the type II seesaw, the proton decay 
rate has been calculated in [ |106|| . In our analysis, we make no such a specific choice, and 
perform detailed analysis in general situations of the minimal SO(10) model by varying the 
above five free parameters. The result for tan/3 = 2.5 is presented in Figure |2]7]. Here the 
distributions of the proton lifetime (log years) for arbitrary choices of the five free parameters 
(normalized by 1) is depicted. We can see that some special sets of the free parameters can 
result the proton lifetime consistent with SuperK results. In that region, cancellation in 
the second line in Eq. (|2.142|) occurs by tuning of the free parameters in the Higgsino mass 
matrix. Note that number of free parameters is not enough to cancel both of the process 
p — > K + T\ (dominant mode) and p — > K + T7^ (sub-dominant mode), and thus the proton 
lifetime has an upper bound in the 50(10) model. For tan = 10, we obtain the same figure 
depicted in Figure [TH] but the lifetime is scaled by roughly (2.5/10) 2 , which is consistent 
with the naive expectation that the lifetime is proportional to 1/ tan 2 /3. Whole region is 
excluded in the case with tan j5 = 10. 

In the case of nondegenerate masses of M^, the parameters increase from three to 5 + 1 
(the last 1 is the ratio of masses). However, the results only slide by the square of this mass 
ratio and do not show the special cancellation. 
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Figure 2.7: The distributions of the proton lifetime (log years) for tan/3 = 2.5 in arbitrary 
five parameter choices (normalized by 1). The green region (right side from 33.3 point) is 
consistent with experiment. 

We have found that for tan /3 = 2.5 some special sets of the parameters predict the proton 
decay rate consistent with the SuperK results, where the cancellation for the dominant modes 
of the proton decay amplitude occurs by tuning of the parameters. Although there exists 
the allowed region, it is very narrow. Our result is consistent with the one in the previous 



work ||106|| for only one specific choice of the Yukawa coupling matrices. It has been found 
that the resultant proton decay rate is proportional to tan 2 (3 as expected and the allowed 
region eventually disappears as tan/3 becomes large, even for tan/3 = 10. 

There are some theoretically possible ways to extend the proton lifetime. One way 
is to adopt a large mass hierarchy between the sfermions and the Wino as can be seen 
in Eq. ( |2.142[) . The proton lifetime is pushed up according to the squared powers of the 
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Figure 2.8: Same as in Figure |2.7| , but for tan (3 = 10. There is very few region (right side 
from 35.3 ) allowed. 



mass hierarchy, and the allowed region becomes wide. How large the hierarchy can be 
depends on the mechanism of the SUSY breaking and its mediation. When we assume the 



minimal supergravity scenario, the cosmologically allowed region ||107|| consistent with the 
recent WMAP satellite data |33| suggests that the masses of sfermion and Wino are not so 
hierarchical and the value we have taken in our analysis seems to be reasonable. Another 
way to evade fast proton decay is to abandon the assumption of Higgsino degeneracy at the 
GUT scale, and to make the mass eigenvalues of the effective colored Higgsino mass matrix 
heavy. We can examine this possibility based on a concrete Higgs sector. However, this 
seems to be a very difficult task even if we introduce a minimal Higgs sector in the minimal 



SO (10) model discussed in ]25] [2(| since there are lots of free parameters in the Higgs sector. 
Furthermore, even in the minimal Higgs sector, there are lots of Higgs multiplets involved 
and the beta function coefficients of the gauge couplings are huge. It seems to be very 
hard to succeed the gauge coupling unification before blowing up of the gauge couplings. 
Therefore, the assumption that all the Higgs multiplets are degenerate at the GUT scale 
would be natural. 

Consequently our results show the typical properties of SO(10) GUT but are not exhaus- 
tive. Also there is possibility to vary GUT phases /?', 7' of ( |2.48|) generically. 
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General Higgs superpotential Also we may consider the more general superpotential 
for completeness [pTfl . 

W = ^mi$ 2 + m 2 AA + ^m 3 H 2 

+ ^m 4 A 2 + ^m 5 E 2 + ^m 6 D 2 

+ Ai$ 3 + A 2 $AA + (A 3 A + A 4 A) H<$> 

5 6 120 

+ E (\ S E 2 + X 9 A 2 + Ai $ 2 + AnA 2 + A i2 A 2 + A 13 # 2 ) 

+ D 2 (XuE + Ai 5 $) 

+ D{\ 16 HA + \ 17 H$ + (A 18 A + A 19 A) (A 20 A + A 21 A) $} . (2.144) 

Here A = 45, A = 126, $ = 210 and E = 54 -plets. For general coupling constants 
Ai, ■ • • ,A 2 i, mi,-- - ,ms, the solutions with higher symmetries are specified by following 
relations. Solutions with higher symmetries are characterized by: 

1. 5*7(5) x U{l) x and (577(5) x U(l)) mppcd symmetry solutions 

E = v R = 0, 



to. (2.145) 



where e = 1 and £ = —1 correspond to the SU(5) x U(l)x symmetric vacua and 
(SU(5) x £/(l)) flippcd symmetric vacua, respectively. For the concrete matter contents, 
see the nextsection. 



2. SU(5) symmetry solutions 
E = 0, 



$! = ^$ 3 , $ 2 = ^$ 3 , At = -^A 2 , V R ^ 0. 

3. G422 = 517(4) x 51/(2) L x 5C/(2) fl symmetry solutions 

$ 2 = $ 3 = A x = A 2 = v R = 0, 
$ x ^ 0, E ^ 0. 

4. G 3221 = SU(3) C x SU(2) L x SU(2) R x C/(1) B _ L symmetry solutions 

$ 3 = Ai = = 0, 
$1 + 0, $ 2 ^0, A 2 ^ 0, £ ^ 0. 

5. G 42 i = 5*7(4) x SU(2) L x C/(l) symmetry solutions 

$ 2 = $ 3 = A 2 = v R = 0, 
$i^0, A x ± 0, £ ^ 0. 



(2.146) 



(2.147) 



(2.148) 



(2.149) 
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6. G 32 u = SU(3) C x SU{2) L x U(1) R x U{1) B -l symmetry solutions 

£ ^ 0(i= 1,2,3), Ai ± (i = 1,2), £7 ^ 0. (2 ' 150) 



The higher symmetry solutions given in Eqs. (|2.145|) -( [2.150|) lead to the crucial consistency 



checks for all results mentioned before. In order to keep the successful gauge coupling 
unification as usual, it is desirable that all Higgs multiplets have masses around the GUT 
scale, but some Higgs fields develop VEVs at the intermediate scale. More Higgs multiplets 
and some parameter tuning in the Higgs sector are necessary to realize such a situation. 

In addition to the issue of the gauge coupling unification, the minimal SO(10) model 
potentially suffers from the problem that the gauge coupling blows up around the GUT 
scale. This is because the model includes many Higgs multiplets of higher dimensional 
representations. 

According to the line of thoughts from (|1.1|) to (|1 . 2|) , it was natural to consider 

L GUT = L'; en + ^ (2.151) 

up to Mp. Here A3 = O(Mp) and gravitation (spacetime structure) appears as a subdomi- 
nant term. However the blow-up before Mp problem shows that such scheme does not exist 
in its naive sense. 

The minimal SO(10) model also is faced on the fast proton decay |29| 



These facts strongly (though not indipensablly) suggest the presence of extra dimensions, 
which not only solves the above problems but also gives new insights for SUSY breaking 
mechanism [|108|| . 

2.9 Solutions to the Problems of the Renormalizable 
SO(10) GUT- Why Extra Dimension ? 



When we accept the superpotential of (|2.103|) , is the breakdown of the gauge coupling uni- 
fication inevitable ? 



2.9.1 Model modifications in 4D 

There are several approaches to this problem making leave the theory in 4D. 

First let us try to remedy the pathologies mentioned in the previous section preserving 
the principles of the renormalizability but discarding minimality of SO (10) GUT, which is 
to add 120. 

The great advantage of minimal SO(10) model was its high predictivity, implying that 
all quark-leptons mass matrices including Dirac and Majorana neutrinos, are completely 
determined. 
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The reason why the gauge coupling unification is broken is as follows. The renormalizable 
SUSY GUT with Higgs fields of high dimensional representation has many Standard Model 
vacua. However such intermediate energy scale is fixed by only single parameter as was 
shown in (^JT§-(gJT§ and also |10| 



= 3(v - l)(v + l)(2v - l)(v* + 5v - 1) 
C126 8w 6 - 27w 5 + 38w 4 - 70w 3 + 87w 2 - 3k; + 3' 1 ' ' 

where v = j^-. So if we add another Higgs, if we retain renormahzabihty, 120 by virtue of 
which can be free. 

C126 



Since 120 has two SM doublets (1, 2, 2) and (15, 2, 2), mass matrices become ||1 10 



M u = c w M w + c { S M 120 + c 126 M 126 

M d = M 10 + M 120 + M 126 

M D = c 10 M 10 + cgJ,M 120 - 3c 126 M 126 (2.153) 

M e = M 10 + cf 2 ] M l20 - 3M 126 

M L = c L M 126 

M R = c R M 126 . 



Here 

Ji) _ (0+) + (0V) (2) (0 + )-3(0V) (3) ~(0-)-3(C) 
120 " -(0-) + (0'-)' 120 ~ -(0-) + (0'-)' 120 ~ -(0-) + (0'-) ' ( } 

where (cj)±) are expectation values of (1, 2, 2) of 120, a.nd(4>'±) are those of (15, 2, 2) of 120. 

In the original model, 126 takes part of Majorana neutrinos, as well as Dirac Fermions 
( p.29| ). In other word, Y\ 2 % is of 0(1) as Yi to recover the wrong SU(5) mass relation ( p,15| ). 



The mass of heavy right handed Majorana neutrino is surely several orders smaller than 
M GUT (we recognised that type II seesaw is subdominant), which means that we are forced 
to have the vev vr of intermediate energy scale. However, we have additionally many pa- 
rameters and can use 126 for determining Mr and Ml independently on the determination 
of Dirac fermion mass matrices. That is Y\ 2 q is free from order one unlike the minimal case 
and vevs are free from having the intermediate energy scales and we can recover the gauge 
coupling unifications. This seems to be fine at least for data fittings of low energy. This 
model has been extensively discussed, especially on suppression of proton decay, by [ [Lll| 

In order that such theory becomes the SM of next generation, we must also study Doublet- 
Triplet (D-T) problem and SUSY breaking mechanism. We will see this point soon later. 

One of the other approaches is to use Split Susy [ |112j| with light gauginos and higgsinos in 



100 TeV range and superheavy squarks and sleptons in energy scale close to GUT. However, 
it is essentially non SUSY and seems to be unnatural. 



The other, for instance, is to adopt non SUSY SO(10) GUT with 45 + 126 Higgs [ |113|1 . 
However, it seems to us these are too conservative to construct new model BSM. 
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2.9.2 Flipped GUT model 

Before we consider flipped SO(10) model we go back to flipped SU(5) model. In SU(5) 
GUT, the SM contents are embedded in different multiplets of 10 and 5. Then hypercharge 
assignment is not unique in these two sets. This may easily understood if we consider two 
series of symmetry breaking (see the results of the previous section) 

50(10) D SU(5) x U(l) v , 3U(5) D SU(3) C x SU(2) L x U(l) z (2.155) 

and 

50(10) D SU(4)p S x SU(2) L x SU(2) R D SU(3) C x SU(2) L x SU(2) R x U(l) x . (2.156) 
In Eq. (|2.155 ), there are two assignments of hypercharge: 

Y - = Z (2.157) 
I = -i(Z + V). (2.158) 

Eq.( |2.157D corresponds to the usual Georgi-Glashow, and Eq. ([2.158|) does to the flipped 
SU(5) [ p.l4j| given from the Georgi-Glashow model by the interchange of SU(2)p doublets 

u c <^d c , e c ^v c (2.159) 

and we obtain harmless relation 

M u = M v (2.160) 

in place of ( |2.15| ) of Georgi-Glashow SU(5) model. Moreover, it is attractive from D-T 
splitting: Higgs superpotential has the form 

W H = 10 x 10 x 5 + TO xlO x 5, (2.161) 

which gives rise to triplet mass 

((1, 1, 0) 10 )(3, 1; 1/3) 10 (3, 1; -l/3) 5 + ((1, 1; 0)^(3, 1; -l/3) w (3, 1; 1/3) F (2.162) 

but has no doublet mass since 5 + 5 has no partner in 10 + 10 (Missing partner mechanism). 

This is a solution to the D-T problem without additional adjoint Higgs. 

However, in flipped SU(5) corresponding to the assigment ( |2.170| ), v c does not belong to 
SU(5) single but to 10-plet, which drives us to unrenormalizable heavy Majorana neutrino 
mass term, 

lOilOjIOtflOtf /M P , (2.163) 

responsible for seesaw mechanism. 

For SO(10) case, SM matter contents are embedded in single 10;. So if we consider 
flipped SO(10), we are forced to enlarge group like 5O(10) x U{l)y C Eq. 
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This implies the addition of matters other than SM matters since 

27 = 16i + 10_ 2 + 1 4 . (2.164) 
There are three kinds of assignments in SO(10) ||115 |: 

16 = d c + e + v + u c + u + d + e c + u c 

10 = D + E c + N c + D c + E + N (2.165) 
1 = S, 



16 = D c + E c + N c + d c + u + d + u c + S 

10 = D + E + N + u c + e + u (2.166) 
1 = e c , 



16 = D c + E + N + u c + u c + u + d+e c + S 

10 = D + E c + N c + d c + e + v (2.167) 
1 = v c . 

These sets are each classified into two different sets and therefore six kinds of classification 
in SU(5). For the first set 

16 = (d c + e + v) + [u c + u + d + e c ) + v c 

10 = (D + E c + N c ) + (D c + E + N) (2.168) 
1 = S 



and 



16 = {u c + e + v) + (d c + u + d + v c ) + e c 

10 = (D + E + N) + (D C + E C + N C ) (2.169) 
1 = 5. 

Flipped SO (10) is given from one of the second set 

16 = {d c + E c + N c ) + {D c + u + d + S) + u c 

10 = (D + e + v) + (u c + E + N) (2.170) 
1 = e c . 



Unfortunately, there is no renormalizable terms making work of Missing Partner mecha- 

(2.171) 



nism unlike flipped SU(5). In this case the counter example is [116] 

W H D T6i16 2 16 2 T6i/M p + 16iT6 2 T6 2 16i/M P 
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W H 3 (1)(10)10 5/Mp + (1)(10)10 5/Mp (2.172) 
which corresponds to ( p,161| ), leading to massive triplets and massless doublets. However, 



we can not incorporate this interaction naively in E 6 model since 

27i27i27 2 27 2 (2.173) 
includes ( |2.172j ) with vev of 16i but simultaneously it also includes 

5i5 2 T0 2 and 5 T 5 2 10 2 , (2.174) 
which give mass to doublet with vev of 16 2 . 



Another approach is to use Dimopoulos-Wilczek mechanism ||117|| . Bertolini et al ||118 
considered flipped SO (10) <S> U(l) model with 45 © 2 x (16 © 16) Higgs fields whose Higgs 
superpotetial is given by 

W H = |Tr45 2 + pijieJSj + ^16,4516, (i,j = 1, 2), (2.175) 

leading directly from 5*0(10) © U(l) to the SM. This sounds good, however, they are forced 
to introduce unrenormalizable Yukawa coupling like 

W Y = YulQ F lQ F lQ H + -^-[Y E 10 F 1 F 1Q H 1Q H + Y d 1Q f 1G f 1G h TQ h ] (2.176) 

to achiev a realistic texture. It seems to spoil the very merits of SO (10) GUT. 

We have discussed D-T splitting in both missing partner mechanism and Dimopoulos- 
Wilczek mechanism. However, there exists no arguments for why this is necessarily the case. 
To assert that D-T splitting is completely solved, we must also explain why \i term is so 
small, which requires a symmetry, most probably R-symmetry in SUSY |30|| . Before we 
address this problem, we argue briefly in the next three subsections on GUT from different 
angles not mentioned so far. 

2.9.3 Perturbative SO(10) GUT 

Our renormalizable model use high dimensional Higgs like 126 and 210, which makes the 
unified gauge coupling blow up after GUT scale, probably before the Planck scale. However, 
GUT scale is O(10 16 )GeV rather near to the reduced Planck scale O(10 18 )GeV, around 
which the renormalization may lose its conventional meaning. Anyhow many physicists 
prefer to adopt low dimensional Higgs fields, which assures the validity of perturbation to 
the Planck scale and is called perturbative SO(10) GUT. Such perturbative GUT leads us 
inevitably to unrenormalizable Yukawa coupling and lose a definite criterion to construct 



Yukawa couplings. As a typical example let us consider the case of Raby's model [|119 



W = W f + W u , (2.177) 
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where 



W f = 16 3 10 H 16 3 + 16 a 10 HX a 

+ Xa ^MxXa + 45 ^ 16 3 + 45 ^ 16a + A16 ^j ' ( 2 - 178 ) 

W v = T6(X 2 N a 16 a + X 3 N 3 16 3 ) + ^(S a N a N a + S 3 N 3 N 3 ) (2.179) 

with a = 1,2. Yukawa coupling unification at Mqvt is realized only for the third generation. 
Their number of parameters is 24, whereas ours is 17 (see subsection ( |2.4.2[ )). Also their 
number counting is quite differnt from ours. Their mass matrices are approximated as real 



but ours are generic. As we pointed out in [fI7 |, complex phases are very important even in 



matching up mixing angles as well as in CP violating processes. Also it is very difficult to 
construct Higgs superpotential generically in the case of perturbative SO(10) GUT case. 

2.9.4 Linear and nonlinear realization 

In the previous subsection we have discussed the larger group E$ and might discuss larger 
group E-j or Eg. The fields transform linearly under the corresponding gauge group, which 
are called linear realization. Let me consider a symmetry breaking from G to H, there 
appear NG bosons. If we consider their superpartners as the SM matters and if we consider 
Lagrangian in terms of only NG bosons, Lagrangian is automatically invariant. Let me 



explain it ||120|| . In this case, the numbers of NG bosons are equal to dim(G/H)=dim G-dim 
H and transforms linearly under h G H transformation but nonlinearly under g G G. Let us 
describe the representatives of G/H as £(7r), parametrizing in terms of ir a as 

£(7r) =<><*(*)//* 9 tt(x) = n a {x)X a , (2.180) 

a&G-H 

where 7r are NG bosons. 

9^) = ^')h(7r,g), h(-K,g)eH (2.181) 

NG field ir(x) is transformed as 

£{n)^£(*')=g£(n)h- 1 (n,g), (2.182) 

under the global g G G transformation. Since this transformation from 7r to tt' is nonlinear 
w.r.t. 7r and the representation of G in terms of tt is called nonlinear realization. The low 
energy effective Lagrangian on G/H is constructed as 

f 2 

C = ^-Tr(d^U j d»U), (2.183) 

where 

U = K T . (2.184) 
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k — 1 


k = 2 


k = 3 


k = 4 


(10, 1/2) 


(10, 9/20) 


(10, 9/22) 


(10, 3/8) 


(16, 5/8) 


(16, 9/16) 


(16, 45/88) 


(16, 15/32) 




(45, 4/5) 


(45, 8/11) 


(45,2/3) 




(54, 1) 


(54, 10/11) 


(54, 5/6) 






(120, 21/22) 


(120, 7/8) 






(144, 85/88) 


(144, 85/96) 








(210, 1) 



Table 2.4: Unitary, potentially massless representations of 50(10) realized at affine levels 
k = 1,2,3,4. Each representation R is listed as (dim R, /ifl) where Iir is its conformal 
dimension. Singlets and conjugate representations are not explicitly written, but understood 
1122 1 . 



Thus even if we consider the same higher gauge group G, we have different low energy 
contents depending on whether we adopt linear or nonlinear realization. In the case we see 
quite different pattern of symmetry breaking i?g or £7 £4 x U(l) 4 from the former case. 
Indeed, Kugo-Yanagida showed that three generations 3 x (10 + 5* + 1) + 5. Unfortunately it 



includes an additional one 5-plet [|121|| . This is quite interesting because it gives new insight 
on the relation between gauge symmetry breking patterns and family symmetry. 

Thus even if we fix gauge group, we have different aspects depending on linear or nonlinear 
representation. 

2.9.5 Constraints from the string theory 

We have discussed GUT so far from the bottom-up approach. It is needless to say that 
top-down approach is also important. 

There are arguments that the string theory does not allow high dimensional Higgs ||122 



This is concerned with Heterotic string model and due to perturbation. Non perturbative F 
theory is out of this constraints, and 126 does not necessarily denied. However, it may be 
useful to consider such counter example. 

As we have shown 126 takes very important roles in renormalizable minimal SO(10) 
GUT model. However, one would require SO(10) at Kac-Moody levels k > 1 (see Table 
2~4j Figure caption is described in terms of string theory terminologies. Affine levels are 



rank of SO(10) in our language. 10, 45, 120, 210, and 126 are rank 1, 2, 3, 4, and 5 
antisymmetric tensor, respectively. 54 is symmetric rank 2 tensor |K) . 
The central charge of SO (10) at level k is given by 

c = Mim{G) = (2.185) 

k + h G (fc + 8) V ; 

for SO(10). Here h G is coxeter. The perturbative heterotic string central charge must be 

c < 22. (2.186) 
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It goes from (|2.185|) and ( p,186| ) that conformal anomaly free condition leads to k < 4. 



Thus if we assume GUT is top-downed from heterotic string [ 123 1 perturbatively, we can not 
produce 126. 

F-theory GUT tries to predict masses of quark-leptons and magnitude of mixing matrices 
in close accord with experiments ||124 . 



There are few papers which have tried to unify bottom-up and top-down approaches 
and these approaches will become more important hereafter. 



rather closely [125 



Then we will go back to the main flow in the next subsection. 

2.9.6 SUSY breaking in 4D 

The tree level potential is given by 



Here 



V 



F F f 



Fl 



dW 



D A = g A Ql(T A r b Q b 



(2.187) 



(2.188) 



dQV ~ a W 
So SUSY is spontaneously broken unless all F a = 0, D A = 0. 

F„ / (called O'Raifeartaigh mechanism) is dominant in many cases and we will con- 
sider mainly this case. D A ^ (called Fayet-Iliopoulos mechanism) is considered only in 
anomalously mediated symmetry breaking in this review. If SUSY is spontaneously broken 
at tree level it must satisfy 



str(M 2 ) = tr(M 2 ) - 2tr(M\ /2 M 1/2 ) + 3tr(Mj 



0. 



(2.189) 



This directly leads us to the result that the scalar superpartners must be lighter than the 



heaviest observed fermion, which contradicts with observation [[72]. So we must consider 
some non-renormalizable process or loop corrections. 

We must separate visible sector from invisible sector where SUSY is broken. There are 
mainly two SUSY breaking mechanisms, gravity mediation and gauge mediation. In either 
case, gravity mediates both sectors. For F-term breaking, gravitino mass is 



"13/2 



(F) 



y/3M P ' 

However, the value (F) is quite different in these two mechanisms. 



(2.190) 



Gravity mediation mechanism The gravity mediation mechanism ||126|| is to consider 
gravitational interaction which communicates SUSY breaking in breaking sector to visible 
sector [ 127|| in 4D. We assume here that SUSY is broken in the hidden sector by the F 
component of a field X. Taking the fact that ultra-violet gravity plays the role into consid- 
eration, we may relax in this case the renormalizability so far assumed. The most general 
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interactions X and the visible sector's fields are [|108 



AC 



d*6 



M 2 



+ 



+ 



Mi 



-XH u H d 



3 x^xq\q 3 + 

V 



Mf 



X^XH u H d + /i.e. 



(2.191) 



d 2 # 



d 2 e 



M P 



■xw?w la + 



+ h.c. 



XQ l H u (u c y + .. 



and we find SUSY breking soft mass, 



m 



soft 



Mp 



(2.192) 



m so f t is of 0(1) TeV and yjF) w O(IO) 11 GeV. This value is very impressive since it 
seems to be common with the breaking energy scales of B-L and Peccei-Quinn symmetries.lt 
should be remarked that b and b' terms give /i and Bfi terms, and that Sj does gaugino 
mass. This model can explain why // and Bfi are suppressed compared with (F). However, 
it simultaneously accompanies flavour problem unless the off diagonal part of |ay| is highly 
suppressed. To assume diagonal is unnatural unless there are flavour symmetry at Planck 
scale. Such symmetry at Mp is unlikely from Black Hole evaporation 
minimal SUGRA ansatz 



108|| . In this sense 



(*?)! 



Z 



(2.193) 



is unnatural. 



Gauge mediation mechanism Another mechanism is to assume that SUSY breaking in 
invisible sector is communicated to the visible sector by heavy chiral fields, called messengers, 
through ordinary SU(3) C x SU{2)l x U{1)y gauge interactions (GMSB). Contrast to gravity 
mediation, it is natural that GMSB is flavour blind at Mp since gauge theory does not 
discriminate flavour. 



Gauge mediation [ 127|| , ||128j| , ||129|| suffers from the anomalously small gaugino masses 



compared to the scalar masses. 

where N m — b — b'. Here b and b' are the beta functions of MSSM and MSSM+messengers, 
respectively. Also it has \x — Bfi problem. 



rn 



mi tan 2 /? - m 2 H 



= -W 2 ^ — ^ ^ (2.195) 

2 IPI tan 2 /3-l V ; 



2Bfi 

2\^\ 2 + m 2 Hu + m 2 



sin2/3 = o| 2 \ 2 ■ ( 2 - 196 ) 
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If we consider /i is generated by the SUSY breaking, 



W = XXH u H d . 

Here X is the backgroud chiral superfield, and X = M + 9 2 F. Then we obtain 

1 F 



H — AM 



-, Bfi = XF « 167r z AM/i = 16vrV > /" 



16tt 2 M„ 



so/i 



47T M mess 47T 



MSSM scalar masses are 

mi = 2A 2 



"3 
47T 



^•) + (£) 2 ^) + (g) 2 ^) 



(2.197) 

(2.198) 
(2.199) 

(2.200) 



If we go beyond 4D and go to extra dimensions, we have two new SUSY breaking mechanism, 
gaugino mediation ||130|| and anomaly mediation (AMSB) ||131|| ||132|| . We will discuss on these 
mechanisms in the next chapters. 

Above the SUSY breaking mechanism, the initial condition of SUSY breaking is also very 
important, which will be discuused in the last subsection of the LHC results. 

2.9.7 No-go theorem in 4D 

There is arguments that it is impossible to construct a GUT in 4D with a finite number of 
multiplets that leads to the MSSM with a residual R symmetry whose no-go theorem is 
not applicable to extra dimensions. This is very important, and let me explain this: SUSY 
invariant action is assumed to be invariant under global U(1)r transformation (for N=l 
supersymmetry as we consider in this review), 

e ->■ e ia e, 9* -> e~ ia e\ (2.201) 

impling that R-charge of 9 and 9^ are 1 and -1, respectively. Chiral superfield is expressed 

$ = <j){y) + V29iP(y) + 99F(y) (2.202) 



as 



with 

•f = x» + i0V0. (2.203) 
Vector superfield is real and its R-charge =0. Vector superfield in Wess-Zumino gauge is 



V = 9^9 A„ + eWd\ + 666* \* + -666*6* D 



(2.204) 



and Ap, X, D have R-charge 0,1,0, repectively. 

Nelson and Seiberg discussed the relation between R symmetry and SUSY breaking [30 . 
They showed under the condition 
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i) Superpotential is generic, and 



ii) low energy theory can be described by a supersymmetric Wess-Zumino model 
that 

a) R symmetry is necessary for SUSY breaking, and 

b) spontaneous R symmetry breaking is sufficient for SUSY breaking. 



Concretely speaking, this will be explained as follows |SD|. Let us consider N superfields 
whose U(1)r charges are -R($) = fa. The R charge of superpotential is 2 and at least one 
superfield must have nonzero R charge, which we take it as N'th field $tv, 0tv 7^ 0, without 
losing generality. If W is R symmetric, then W can be described as 



W($i) = (2.205) 



where 



**=j!fc> « = 1.-^- 1. (2-206) 



•IT,' I — 'I'v' V 

dW </w N 



(2.207) 



^2/<t> N BW 

N 



Now if U(1)r are spontaneously broken, ($jv) 7^ 0, ( [2.207D leads to 



~ dW 

W = and = 0. (2.208) 

Generically we can not satisfy the above N constraints by iV — 1 fields. 

Thus if we have no U(l) symmetry we have appropriate SUSY vacuum, that is, U(l) 
symmetry is necessary for SUSY breaking (condition (a)). 

If there is U(l) symmetry and it is spontaneously broken, SUSY is automatically broken 
(condition (b)). 

So the problem is how to impose U(1)r symmetry in superpotential of GUT. 

Reflecting these situations, Fallbacher et al. ]3T] concluded that no MSSM model with 
either a Zf f>3 or U(1).r symmetry can be completed by a four dimensional GUT in the 
ultraviolet. The essential point is explained for SU(5) GUT as follows. SU(5) x is 
broken to the SMxZ^ by the vev of the SM singlet of 24. 24 has zero R-charge since 
is unbroken, and 



24 = (8, 1) © (1, 3) © (1, 1) © (3, 2)_ 5/6 © (3, 2) B/fl . (2.209) 

Here ((1, 1) ) 7^ 0, and (3, 2)_ 5 ^ 6 and (3, 2) 5 y 6 get absorved to the longitudinal part of gauge 
bosons. The remaining (8, 1) and (1,3) must be massive and therefore require mass term 
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m24 x 24. However, it is prohibited because 24 have R-charge but superpotential must 
be 2 R-charge. this is the case for more general mutiplet and more general gauge group 
including SO(10), The details should be referred to [|3T|. On the other hand in the case of 
Pati-Salam case, PS group to the SM need to reduce rank by one, which is done by (4,1,2) 
and break B-L quantum number and there give rise to no problem. Therefore, the minimum 
group subject to no-go theorem is SU(5). 

Of course there are a loophole of this no-go theorem. For instance it is for meta-stable 
supersymmetry breaking vacuum, where U(1)r is broken explicitly |133 |. That is, let us 
consider 

W = -Jfe$i +m$ 2 $ 3 + |$i$3, (2.210) 

which is U(1)r symmetric with R-charge, = R$ 2 = 2, R§ 3 = 0. They introduced a 
broken term, 

AW = ^em<5>l (2.211) 

where e is a small dimensionless parameter. However, this time we must explain why e is so 
highly tuned to satisfy longevity of metastable state $i = $ 2 = $3 = and we do not adopt 
this scenario. 

On the other hand, no-go theorem can not be applied in an extra dimensions, where new 
ways of GUT symmetry breaking mechanisms appear [|134|| [ 135 ] [32| . This is one of very 
strong motivations for us to proceed to extra dimension. []. 



We may consider ( |2.201| ) from string theory. In string theory [|137|| , it has originally global 



space-time SO(10) symmetry and is broken to 50(4) x SO (6) in 4D. This SO(6) is isomorphic 
to SU(4). The spinor in ten space-time dimensions has 16 l + 16r components. (Do not 
confuse with flabour group so far discussed.) In the splitting from 10 to (4+6) dimensions, 
this spinor is divided into four 4-component spinor, 9a\0^ (a = 1,2), i = 1,2,3,4 So 
there is S77(4)# transformation 

0/(0 = U i d (j) m (2.212) 

Taking all the arguments on the problems in the minimal SO(10) into considerations, 
we will mainly study the possibility of extending GUT to extra dimensions in the susequent 
sections. 

Going beyond 4D to 5D or more extra dimensions, we achieve much more merits than 
demerits as we will explain. All previous correct predictions are remained valid, whereas 
the small mismatch comes from the tight mass relation of ( f2.29| ) will be improved. Gauge 
coupling unification is recovered. Fast proton decay is evaded. 

Though we have discussed bottom-up approach, it is useful or indispensable to consider 
top-down approach also |6]| |[L22|| . In this the scale of the bottom-up approach is close 

to Planck scale, it is expected to coincide with top-down approaches from string world J33 



The mutual coincidence and compatibility give also some informations to both approaches. 



7 The roles of more general no-go theorems in GUT were discussed in 136 
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Comment against the Eingorn-Zhuk's arguments Before discussing SUSY GUT in 
5D, we discuss the Eingorn-Zhuk's argument | 138 | since many physicists consider it as the 
no-go theorem against the presence of the extra dimension. 

They considered the reduced action in D+l dimensions under the assumption that 

//;, >'/,, (hj = 0,1, ...,£>) (2.213) 

at spatial infinity, 

S = —Et + Mijj + S 3 (r 3 ) + S 4 (x 4 ) + .... + S D (x D ). (2.214) 

Here E and M are conserved energy and angular momentum of the system, respectively. 
However, this action is consistent only if S is independent on the coordinates of extra di- 
mensions, x 4 , ...,x . In that case only, Mip term can be consistent and conserved quantity 
since angular momentum in general defined by 

= ^Gt/V - pV). (2.215) 

In their model, spatial isotropy is assumed only in usual three dimensions and therefore test 
particles have zero momentum 

= (2.216) 

in extra dimensions, which means 

d IM S = 0f<3rfi, = 4,..,D. (2.217) 
Under these assumptions they obtained 

n, s f\f r, E2 \ 1^ 2 2 2(£> - 1) „ \ 1 / „ r2 Dm 2 c 2 r 2 q \ 1 1/2 

(2.218) 

where r g = with the usual Newton constant G and source mass M. For the Mercury 
perihelion precession is given by [|139 



dr 



DTTm 2 c 2 r 

** = W^W*' (2 ' 219) 

which is consistent with observation if D = 3, irrelevantly to the size of the extra dimensions. 
It follows from these arguments that the extra dimensions are excluded if 

• Space-time is asymptotically flat, 
and 

• Fields other than gravitation does not enter into extra dimensions. 
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Of course, GUT models in extra dimensions which will be discussed in the next chapter are 
free from these constraints. 

We consider SO(10) is realized in 5D and the symmetry in 4D is replaced by the Pati- 
Salam invariance. As the result of this change, the mass matrix of heavy right-handed 
neutrino which was represented by M 12 6 ( P-29Q , one of the partners of mass matrices of 
charged fermions and Dirac neutrinos ( |2.29| ), is replaced by the independent mass matrix 
( p. 20 ). For that reason, mismatch of the data is resolved. Also we can eliminate the Higgs 



component which was involved in 126 and induced fast proton decay. However, this is not the 
all reasons why we must go beyond 4D. We can recover the gauge coupling unification spoiled 
by the intermediate energy scales. Furthermore, it is free from the internal inconsistency like 
no-go theorem on the SUSY breaking mechanism. 
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Chapter 3 

Orbifold SO(10) GUT in Five 
Dimensions 



In this chapter we exploit the orbifold SO(10) GUT. By the compactification on the orbifold 
32], N — 1 SUSY of the five dimensional theory, corresponding to iV = 2 SUSY in the four 
dimensional point of view, is broken down to four dimensional N = 1 SUSY. Supersymmetric 
Lagrangian of this system can be described in terms of the superfield formalism of four 
dimensional N 



1 SUSY theories 140, 141, 142 



There are lots of possibilities to construct 
such a model, where some fields reside in the bulk and some reside on branes. Bulk itself 
has two possibility that it is flat or warped. We first discuss the case of warped model and, 
subsequently main flat case, orbifold GUT. 



3.1 Warped SO(10) GUT 



We embed the minimal SO(10) model [ 143 into a warped extra dimension model [F4"3|. In 



this scenario, the warped metric gives rise to an effective cutoff in 4-dimensional effective 
theory, which is warped down to a low scale from the fundamental mass scale of the original 
model (a higher dimensional Planck scale). We choose appropriate model parameters so as to 
realize the effective cutoff scale as GUT scale. Furthermore, in the context of a warped extra 
dimension we can propose a simple setup that naturally generates right-handed neutrino 
masses at intermediate scale even with Higgs field VEVs at the GUT scale. Thus, the gauge 
coupling unification remains as usual in the MSSM. The fifth dimension is compactified on 
the orbifold S 1 /Z2 with two branes, ultraviolet (UV) and infrared (IR) branes, sitting on 
each orbifold fixed point. With an appropriate tuning for cosmo logical constants in the bulk 



and on the branes, we obtain the warped metric | 144 | 



ds 2 = e- 2kr ^r lpLV dx^dx v - r 2 c dy 2 , (3.1) 

for — 7r < y < 7r, where k is the AdS curvature, and r c and y are the radius and the angle of 
S 1 , respectively. 

The most important feature of the warped extra dimension model is that the mass scale of 
the IR brane is warped down to a low scale by the warp factor u = e~~ krc7T in four dimensional 
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effective theory. For simplicity, we take the cutoff of the original five dimensional theory and 
the AdS curvature as M 5 ~ k ~ Mp = 2.4 x 10 18 GeV, the four dimensional (reduced) Planck 
mass, and so we obtain the effective cutoff scale as Ajp = uMp in effective four dimensional 
theory. Now let us take the warp factor so as for the GUT scale to be the effective cutoff 
scale M GU t = A IR = uM P , namely uj ~ 0.01 ||145|| . As a result, we can realize, as four 
dimensional effective theory, the minimal SUSY SO(10) model with the effective cutoff at 
the GUT scale. 

Before going to a concrete setup of the minimal SO(10) model in the warped extra 
dimension, let us see Lagrangian for the hypermultiplet in the bulk, 



d 2 ee -3kr c \y\ H c 



\ - c) kr c e(y) - j= 2 



H + h.c. 



(3.2) 



where c is a dimensionless parameter, e(y) = y/\y\ is the step function, H, H c is the hyper- 
multiplet charged under some gauge group, and 



V 

X 



-Qa^QA^ - %9'9X 1 + i9 z 9X 1 + -9'9 Z D 

-^=(E + iA 5 ) + V29X 2 + 9 2 F 
v2 



(3.3) 



are the vector multiplet and the adjoint chiral multiplets, which form an iV = 2 SUSY gauge 
multiplet. Z2 parity for H and V is assigned as even, while odd for H c and x- 

When the gauge symmetry is broken down, it is generally possible that the adjoint chiral 
multiplet develops its VEV ||146|| . Since its Z 2 parity is odd, the VEV has to take the form, 



(E) = 2akr c e(y) 



(3.4) 



where the VEV has been parameterized by a parameter a. In this case, the zero mode wave 
function of H satisfies the following equation of motion: 



H = 0, 



(3.5) 



which yields 



H 



D (3/2-c+a)kr c 



N 



(3.6) 



where h{x^) is the chiral multiplet in four dimensions. Here, N is a normalization constant 
by which the kinetic term is canonically normalized, 



1 



;i - 2c + 2a)k 



]\f e {l-2c+2a)kr c TT _ I ' 



(3.7) 
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Hence, at y — tt, the wave function becomes 

H(y = tt) ~ v/(l-2c + 2a)A; w -1 /i(x^) (3.8) 

if e ( 1 /2-c+a)fcr c 7r > 1; while 

= tt) ~ ^-(1 -2c + 2a)k w -l e (l/2-c+a)fcr c 7T ^) (3 9) 

for e (V2-c+a)fcr c 7r ^ ^ 



Lagrangian for a chiral multiplets on the IR brane is given by 



C IR = / d A 9 uj ] uj $ f $ + 



J d 2 9 u 3 W{§) + h.c. 



(3.10) 



where we have omitted the gauge interaction part for simplicity. If it is allowed by the gauge 
invariance, we can write the interaction term between fields in the bulk and on the IR brane, 

C mt = [ d 2 6u 3 ^=&H{y = tt) + h.c. , (3.11) 
J V M5 

where Y is a Yukawa coupling constant, and M 5 is the five dimensional Planck mass (we 
take M 5 ~ Mp as mentioned above, for simplicity). Rescaling the brane field $ — > <§>/uj to 
get the canonically normalized kinetic term and substituting the zero-mode wave function 
of the bulk fields, we obtain Yukawa coupling constant in effective four dimensional theory 
as 

Y AD ~Y (3.12) 

for e (l/2-*+-a)*rcir ^ ^ wMq 

F 4D ~yxe (1/2 - c+a)tw «r, (3.13) 

for e (V2-c+a)*r c 7r ^ L In the latter 

case, we obtain a suppression factor since H is localized 

around the UV brane. 

Now we give a simple setup of the minimal SO(10) model in the warped extra dimension. 
We put all 16 matter multiplets on the IR (y = tt) brane, while the Higgs multiplets 10 
and 126 are assumed to live in the bulk. In Eq. (|3.11|) , replacing the brane field into the 
matter multiplets and the bulk field into the Higgs multiplets, we obtain Yukawa couplings 
in the minimal SO(10) model. The Lagrangian for the bulk Higgs multiplets are given in 
the same form as Eq. Q3.2p, where \ ls the SO(10) adjoint chiral multiplet, 45. As discussed 
above, since the SO(10) gauge group is broken down to the SM one, some components in \ 
which is singlet under the SM gauge group can in general develop VEVs. Here we consider a 
possibility that the U(l)x component in the adjoint x = 45 under the decomposition SO(10) 
D SU(5) x U(l) x has a non-zero VEVf], 

45 = In © 10 +4 © 10_4 © 24n . 



1 Since \ has an odd Z2 parity, its non-zero VEV leads to the Fayet-Iliopoulos D-terms localized on 



both the UV and IR branes [ 147 1 , which should be canceled to preserve SUSY. For this purpose, we have to 
introduce new fields on both branes by which the D-terms are compensated. If such fields are in the same 
representations as matters or Higgs fields like 16 or 126, we would need to impose some global symmetry 
to distinguish them. 
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The 10 Higgs multiplet and the 126 Higgs multiplet are decomposed under SU(5) x U(l)x 
as 

10 = 5+2 © 5_ 2 , 
126 = l+io © 5 +2 ©T0 +6 ©15_ 6 ©45_2© 50+2 . 

In this decomposition, the coupling between a bulk Higgs multiplet and the U(l)x component 
in x is proportional to U(l)x charge, 

dm D^J d 2 6io 3 Q x (i: x }H c H + h.c. , (3.14) 

and thus each component effectively obtains the different bulk mass term, 

Q-cW + iQ x (£ x }, (3.15) 

where Qx is the U(l)x charge of corresponding Higgs multiplet, and Y*x is the scalar compo- 
nent of the U(l)x gauge multiplet (lo)- Now we obtain different configurations of the wave 
functions for these Higgs multiplets. Since the l+i Higgs has a large U(l)x charge relative 
to other Higgs multiplets, we can choose parameters c and (£x) so that Higgs doublets are 
mostly localized around the IR brane while the l+io Higgs is localized around the UV brane. 
For example, the parameter choice, c = —7/2 for both 10 and 126 Higgs multiplets and 
(Ex) = — kr c , can realize this situation. 

Using the decomposition of matter multiplets, 

16* = rl 5 © 5^010!.! , 

the Yukawa couplings between matters and the 126 Higgs multiplet on the IR brane are 
decomposed into 

W Yl26 = 1^5 +2 10L 1 10i 1 + l?45_ a 5 i + al0i 1 
+ F^5+ 2 li 5 5i3 + n ii 45-25t 3 10i 1 

+ + ^l+ 10 li 5 li 5 • (3.16) 

Here, all the Yukawa couplings coincide with the original Yukawa coupling Y" 12 6 up to ap- 
propriate CG coefficients. As discussed above, the l+io Higgs multiplet giving masses for 
right-handed neutrinos is localized around the UV brane and, therefore, we obtain a sup- 
pression factor as in Eq. ( |3.13|) for the effective Yukawa coupling between the Higgs and 
right-handed neutrinos. In effective four dimensional description, the GUT mass matrix 
relation is partly broken down, and the last term in Eq. fl2.28p is replaced into 

Yii 6 v R ^YH 6 (ev R ), (3.17) 

where e denotes the suppression factor. By choosing appropriate parameters so as to give 
e ~ 10~ 2 , we can take vr ~ M GU t and keep the successful gauge coupling unification 
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in the MSSM. In fact, the above parameter set, c = —7/2 and (Ex) = —kr c , leads to 
e = to = Mqut/Mp ~ 10 -2 . The other Higgs multiplets are localized around the IR brane, 
so that there is no suppression factor for other effective Yukawa couplings. 

In our setup, all the matters reside on the brane while the Higgs multiplets reside in the 
bulk. This setup shares the same advantage as the so-called orbifold GUT |H| |148| , |149| . 



We can assign even Z 2 parity for MSSM doublet Higgs superfields while odd for triplet 
Higgs superfields, as a result, the proton decay process through dimension five operators are 
forbidden. 

In order to solve these problems, we have considered the minimal SO (10) model in the 
warped extra dimension. As a simple setup, we have assumed that matter multiplets reside 
on the IR brane while the Higgs multiplets reside in the bulk. The warped geometry leads 
to a low scale effective cutoff in effective four dimensional theory, and we fix it at the GUT 
scale. Therefore, the four dimensional minimal SO(10) model is realized as the effective 
theory with the GUT scale cutoff. 

After the GUT symmetry breaking, the adjoint scalar in the gauge multiplet in five 
dimensional SUSY can generally develop a VEV, which plays a role of bulk mass for the 
bulk Higgs multiplets. This bulk mass is proportional to the charge of each Higgs multiplets 
and cause the difference between wave functions of each Higgs multiplet. We have found 
the possibility that the singlet Higgs which provides right-handed neutrino with masses is 
localized around the UV brane and the geometrical suppression factor emerges in Yukawa 
couplings of the right-handed neutrinos. As a result, we can set the mass scale of the right- 
handed neutrinos at the intermediate scale, nevertheless the singlet Higgs VEV is around 
the GUT scale. All Higgs multiplets naturally have masses around GUT scale and the gauge 
coupling unification in the MSSM remains unchanged. 

We give some additional comments. One can easily extend our setup to put some of 
matter multiplets in the bulk |]15U| , |151| , |152| j. In this case, we may explain the fermion mass 
hierarchy in terms of the different overlapping of fermion wave functions between different 
generations. In this section, we have assumed that GUT gauge symmetry is successfully 
broken down to the SM one. There are several possibilities for GUT symmetry breaking. It 
is easy to introduce appropriate Higgs multiplets and superpotential so as to break the GUT 
symmetry on a brane as in four dimensional SO(10) models. We also be able to introduce 
an appropriate boundary conditions for bulk gauge multiplets to (explicitly) break the GUT 
symmetry to a subgroup with rank five in total, as in the orbifold GUTs. 



3.2 Model Setup of SO(10) GUT in 5D 

From this section we will realize the new model compatible with no-go theorem discussed in 
the last part of previous chapter. That is, we will discuss SO (10) orbifold GUT proposed 
by Ref. P4"| . From the view of string theory, we should consider ten dimensional world. 
However, we will discuss essentially bottom up approach and incorporate the indispensable 
effect of extra dimension as minimally as possible. 

The model is described in 5D and the fifth dimension is compactified on the orbifold 
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S 1 /Z2 x Z' 2 . A circle S 1 with radius R is divided by a Z2 orbifold transformation y — > —y 
(y is the fifth dimensional coordinate < y < 2ixR) and this segment is further divided 
by a Z 2 transformation y' — > —y' with y' = y + irR/2. There are two inequivalent orbifold 
fixed points at y = and y = nR/2. Under this orbifold compactification, a general bulk 
wave function is classified with respect to its parities, P — ± and P' — ±, under Z 2 and Z' 2l 
respectively. 

Assigning the parity (P,P') the bulk SO(10) gauge multiplet as listed in Table I, only 
the PS gauge multiplet has zero-mode and the bulk 5D N=l SUSY SO(10) gauge symmetry 
is broken to 4D N=l SUSY PS gauge symmetry. Since all vector multiplets has wave 
functions on the brane at y = 0, SO(10) gauge symmetry is respected there, while only the 
PS symmetry is on the brane at y = 7tR/2 (PS brane). 



(P, p>) 


bulk field 


mass 


(+,+) 


U(15,l,l), U(l,3,l), U(l,l,3) 


2n 
R 


(+,-) 


U(6,2,2) 


(2n+l) 
R 


(-,+) 


$(6,2,2) 


(2n+l) 
R 


(-, -) 


$(15,1,1), $(1,3,1), $(1,1,3) 


(2n+2) 
R 



Table 3.1: (P, P') assignment and masses (n > 0) of fields in the bulk SO(10) gauge 
multiplet (V, $) under the PS gauge group. V and $ are the vector multiplet and adjoint 
chiral multiplet in terms of 4D N=l SUSY theory. 

Boundary condition in brains are very elegant introduction to GUT model, whose original 
idea was introduced to circumvent the mirror particles. These particles appear when we 
incorporate all matters of three generations preserving chiral symmetry. 

We place the all matter and Higgs multiplets on the PS brane, where only the PS sym- 
metry is manifest so that the particle contents are in the representation under the PS gauge 
symmetry, not necessary to be in SO(10) representation. For a different setup, see The 
matter and Higgs in our model is listed in Table ft.lj For later conveniences, let us introduce 



"Such orbifold condition is also geometrically derived from the singularity free condition 153 
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the following notations: 



H, = (1,2,2V, #( = (1,2,2)'^, 
H 6 = (6,1,1)*, £ib = (15,1,1)h, 
# £ = (4,2,1)*, SZ=J4,2,1)h, 
h, 0= (4,1,2k 



(4,1,2 



(3.18) 



under the PS group ST/(4) C x SU(2) L x S'£/(2) fl . 

Superpotential relevant to fermion masses is given byQ 



Y 



Y?F U F<^H X + ^F U F% (H[H 15 ) 



+ 



R rpc rpc 



(3.19) 



Here the notations are as follows: M5 is the 5D Planck scale. Fn and F^ are matter 
multiplets of i-th generation in (4,2,1) and (4,1,2) representations, respectively. Hi = 
(1,2,2), H[ = (1,2,2)', H 15 =J15, H 6 = (6, 1, 1) H (See (pi])), = (4, 1, 2) H , 

<P = (4, 1, 2) H , H l = (4, 2, 1)*, H L = (4, 2, 1)* are Higgs multiplets. 

The product, H[Hi 5 , effectively works as (15, 2, 2)h, while <p<p and H^Hl effectively work 
as (10, 1, 3) and (10, 3, 1), respectively, and are responsible for the left- and the right-handed 
Majorana neutrino masses. Thus Wy inherits the essential part of the minimal SO(10) 
GUT model. Providing VEVs for appropriate Higgs multiplets, fermion mass matrices are 
obtained. 





brane at y — nR/2 


Matter Multiplets 


i/) i = F Li ®F% i (z = l,2,3) 


Higgs Multiplets 


(1,2,2V, (1,2,2)^, (15,1,1V, (6,1,1V 
(4,1,2V, (4,1,2V, (4,2,1V, (4, 2, IV 



Table 3.2: Particle contents on the PS brane. F14 and F^ { are matter multiplets of i-th 
generation in (4, 2, 1) and (4, 1, 2) representations, respectively. 

3 For simplicity, we have introduced only minimal terms necessary for reproducing observed fermion mass 
matrices. 



60 



M u = c w M 1)2>2 + ci 5 M 15 , 2 ,2 , 

M d = M 1A2 + M 15i2j2 , 

M D = c 10 M li2)2 - 3ci5M 15 , 2 , 2 , (3.20) 

M e = M li2>2 - 3M 15 . 2)2 , 

M L = c L Afio, 3 ,i , 

M R = C R M 10A ,3 . 

Two important remarks are in order. 
1- Mi5,2,2 is, m general, not symmetric unlike Mx 2 6- However, we imposed the L-R symmetry 
4,1,2 -H- 4,2,1, which implies that bothy Mi j2>2 and Mis i2>2 matrices are symmetric and 
mass structure of charged fermions and Dirac neutrino is same as that in SO (10). 
2. Mi and Mr are independent on those of the charged fermions and the Dirac neutrino 
unlike the SO(10) case ( |2.29| ). So the precise data fitting becomes possible without changing 
Y v . This is very important especially for LFV and leptogenesis as will be discussed in the 
subsequent sections. 

Hq is necessary to make the color triplet heavy. However, there arises no D-T problem 
since they are not involved in the same multiplet. There are sufficient numbers of free 
parameters to fit all the observed fermion masses and mixing angles. 

We introduced Higgs superpotential invariant under the PS symmetry |34]] such as Q 

W PS = ^Hf + ^-H[ 2 + miB tr [H 2 15 ] + m 4 (Hj~H L + #) 
+ {H]jp + HE<j>) (Aiffi + KH[) + Ais (# + HlH L ) H 15 
+ \ ti [H\^ + \ & {h 2 l + TT1 2 + ct> 2 + f^H & . (3.21) 



Parameterizing (i?is) = ^%diag(l, 1, 1, — 3), SUSY vacuum conditions from Eq. Q3.21 
and the D-terms are satisfied by solutions, 



8m4 / A . 

wis = ttt— m 4 , ((p) = ((p) = a — y- mis - t— m 4 = wps (3.22) 



3Ai5 y 3A 15 \ A15 

and others are zero, by which the PS gauge symmetry is broken down to the SM gauge 
symmetry. We choose the parameters so as to be vi§ ~ ((f)) = ((f)). Note that the last term 
in Eq. ( ^.21[ ) is necessary to make all color triplets in and heavy. 



Weak Higgs doublet mass matrix is given by 



nil 

H u H[, H L )\ 0_ m' L X'M | [ H[ | . (3.23) 
Ai(0) \'M 










(H H )\ 




\ttJ 



4 



It is possible to co nsid er a different superpotential by introducing a singlet chiral superfield, which will 



be discussed in section 3.4 
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In order to realize the MSSM at low energy, only one pair of Higgs doublets out of the above 
tree pairs should be light, while others have mass of the PS symmetry breaking scale. This 
doublet-doublet Higgs mass splitting requires the fine tuning of parameters to satisfy 

det M = m\rr{ 1 m± — (miAf + m'iXl)vp s = 0. (3.24) 

Suppose a Higgs superpotential which provides the same VEVs (fps) for Higgs multiplets 
to break the PS symmetry and leaves only the particle contents of the minimal supersymmet- 
ric Standard Model (MSSM) at low energies. In Ref. pll , assuming M c = t> PS and imposing 



the left-right symmetry, the gauge coupling unification was examined. Analyzing the gauge 
coupling running in the MSSM, M c = t> PS is fixed as the scale where the SU(2)^ and SU(2)# 
gauge couplings coincide with each other, which is found to be M c = vp$ = 1.2 x 10 16 GeV. 
For the scale \i > M c = vps, we have only two independent gauge couplings, SU(4) C and 
SU(2)l (or SU(2)ij) gauge couplings. After taking KK mode contributions into account, it 
was shown that the gauge coupling is successfully unified at M GU t = 4.6 x 10 17 GeV (see 



Figure 1 from Ref. ||34|| ). We assume that a more fundamental SO (10) GUT theory takes 
place at Mgut-, and it would be natural to assume Mqut ~ M$. In fact, the relation be- 
tween 4D and 5D Planck scales, M|/M c ~ Mp (M P = 2.44 x 10 18 GeV is the reduced Planck 
scale), supports this assumption with M c = 1.2 x 10 16 GeV. When we abandon the left-right 
symmetry, there is more freedom for the gauge coupling unification with two independent 
parameters vps and M c . 

Our model gives a large vps relative to other 5D orbifold SO(10) models |33|]. The high 



value of vps or M c is advantageous for dangerous proton decay due to dimension six opera- 
tors. From Eq. ( p.!9[ ), the right-handed neutrino mass scale is given by Mp ~ Ypvp s /M^ ~ 



Ypvpg/McuT- The scale Mp = O(10 14 GeV) preferable for the seesaw mechanism can be 
obtained by a mild tuning of the Yukawa coupling Y R ~ 0.1. 

In the next section, we show that this model is well fitted to the smooth hybrid inflation 
model and the model predictions are compatible with the cosmological observations like 
the power spectrum of the curvature perturbations, the scalar spectral index, the ratio of 
scalar-to-tensor fluctuations and so on. 

Proton decay Dimension five operator in nucleon decay does not appear since H e does 
not couple with matters as in ( |3.19| ) . Dimension six operator is also suppressed because Mx 
is increased in 5D. 



3.3 Gauge Coupling Unification in 5D 

In the orbifold GUT model, we assume that the GUT model takes place at some high 
energy beyond the compactification scale. For the theoretical consistency of the model, the 
gauge coupling unification should be realized at some scale after taking into account the 
contributions of Kaluza-Klein modes to the gauge coupling running [|154|1 ||155|| . 



In our setup, the evolution of gauge coupling has three stages, G321 (SM+MSSM), G422 
(the PS) and M c = 1/R. For simplicity, we assume vps = M c in our analysis. Furthermore, 
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since we have imposed the left-right symmetry, SU(2)l and SU(2)p gauge couplings must 
coincide with each other at the scale \l — Vps- As a consequence, the PS scale is fixed from 
the gauge coupling running in the MSSM stage. 
In the G321 stage, we have 

1 1 1 , fM\ , , 

-bM — ; (2 = 3,2.1), (3.25) 



ai(n) OLi{M) 2ix \ /i J 
were 6jS are 

63 = -7, b 2 = -19/6, 61 = 41/10 (3.26) 

for Mz < /i < Msusy and 

h = -3, 62 = 1, h = 33/5 (3.27) 

for M S usy < H < M c = up 5 . At the PS scale, the matching condition holds 

a^(M c ) = al\M c ) 
a^(M c ) = a~l(M c ) 

a^(M c ) = [2a^(M c ) + 3a-^M c )}/5 (3.28) 
In the PS stage fi > M c , the threshold corrections Aj due to KK mode in the bulk are added, 

+ — bM ( — - ) + Aj. (2 = 4, 2 L , 2 R ) (3.29) 



ai(fji) ai(M c ) 2ix \ fi 
The beta functions from the matter and Higgs multiplets on the PS brane are 

W = 3, b 2L = b 2R = 6. (3.30) 
KK mode contributions are given by 



Z7T L — ' U 



n=0 



+ T^f* Y OQi - (2n + l)M c )ln (2n + 1)Mc (3.31) 

n=0 ^ 



with 



6f en = (-8,-4,-4), 

b odd = (_ 8j _ 12j _12) (3.32) 

under 6*422- Fig. |3.1| shows the gauge coupling unification for the left-right symmetric case. 
As the characteristic property of gaugino mediation, KK modes affect on the asymptotic 
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Figure 3.1: Gauge coupling unification in the left-right symmetric case, taken from Ref. [46 . 
Each line from top to bottom corresponds to gs, g 2 and g± for ji < M c = vps, while g% = 
and g 2 = g 2 p for \i > M c = v PS . 



free direction. The PS (compactification) scale, M c , is determined from the gauge coupling 
running in the MSSM stage by imposing the matching condition, (M c ) = a^ R (M c ) = 
(5a^\M e ) -2a3 1 (M c ))/3, and we find 



Vps = M c = 1.2 x 10 lb GeV. (3.33) 

for the inputs, («i(M z ), a 2 (M z ), a 3 {M z )) = (0.01695, 0.03382, 0.1176) and M SUS y = 1 TeV. 
For the scale fi > M c , there are only two independent gauge couplings a 4 and a 2 = a 2 R, and 
so the gauge coupling unification is easily realized. We find the unification scale as 

M GUT = 4.6 x 10 17 GeV. (3.34) 

As mentioned before, we assume that a more fundamental SO(10) GUT theory takes place 
at Mquti and it would be natural to assume Mqut ~ M 5 . In fact, the relation between 
4D and 5D Planck scales, M|/M c ~ Mj, (M P = 2.4 x 10 18 GeV is the reduced Planck 
scale), supports this assumption with M c = 1.2 x 10 16 GeV. When we abandon the left-right 
symmetry, there is more freedom for the gauge coupling unification with two independent 
parameters vps and M c . We will consider the gauge coupling unification for 

M c > vps (3.35) 

in section |3.5| . ( |3.35|) is necessary for neutralino LSP in the gaugino mediation. 
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3.4 Confrontation with Cosmology— Smooth Hybrid In- 
flation 



Inflation model is already accepted as the new Standard Cosmology since it solves the short- 
comings of old Standard Big-Bang Model |B6 . Inflation is supposed to have occurred at 
GUT scale and must have strong relationships with GUT. After old inflation | |156| |, there 
appeared many inflation models: slow roll inflation |95j and K-inflation ||157|| are driven by 
the potential and kinetic enery of scalar fields, respectively. Also there are models due to 
modified gravity ||158|| and due to Bose-Einstein condensatation [ 159|| . Inflation model has 
the observational supports from WMAP. Inflation is considered to have occurred at GUT 
scale and its scenario crucially depends on GUT model. Also WMAP data give them the 
severe constraints mainly due to baryon acaustic oscillation (BAO) [|160|1 ||161J| . 

The inflation have three main observables; power spectrum n s , tensor-to-scalar ratio, and 
(non) Gaussianity J38|. Let me introduce them. The spatial elements of the metric tensor 
g^ v are divided into the trace part 1Z and the traceless parts (e' 1 ).^ 



'.i 



9ij = a\t)e™ (e h ) t 
with det(e h ).. = 1. The power spectrum of 1Z is defined by 

and the spectral index n s is defined as 

P n {k) oc k n °-\ 



(3.36) 



(3.37) 



(3.38) 



The current data is n, = 0.96 ± 0.01 



o CL). The tensor-to-scalar ratio r is defined by 
P h (ko) 



Pn(k )' 



(3.39) 



Here Ph{k) is 



(3.40) 



(h ijM ti*> k ) = (27r) 3 6(k + k')P h (k) 

and ko = 0.002Mpc _1 is some reference wavenumber at which r is defined. The current upper 
bound is r < 0.24(95% CL). The three-point function 



(n^TZ^n^) = (2n) 3 5 D (k 1 + k 2 + k 3 )B n {k l7 k 2 , k 3 ) (3.41) 
vanishes for Gaussian fluctuations. Single field inflation is nearly Gaussian, that is, 

6, Bft(fci, k 2 , k 3 ) 

5 JNL Pn(h)Pn(k2) + (2perm) 1 ' ' 

almost vanishes, and non-Gaussian signal indicates the exclusion of it. The current bound is 
f NL = 32 ± 21(68% CL) and n s = 0.96 ± 0.01 (68% CL) g§. The Planck data are expected 
to be open in 2013 p2| . 
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There are some deficits for a single-field inflation model. It gives rise to fine tuning 
problem. For single field inflation f^L is related with n s as 



f NL = ^(l-n s ). (3.43) 

Hence if J'nl > 1 will be observed, it rules out single field inflation. Here we restrict our 
arguments on hybrid inflation [4(| since single field inflation suffers from several tensions 



from observation. The most naive superpotential may be 

W = kS(—M 2 + 00), (3.44) 

where and is a pair of the SM singlet. However, this model is flat along the inflation 
valley (0 = = 0) and must include one-loop correction of the Coleman- Weinber type. We 
want to corelate with high dimensional GUT field. However, and develop vevs during 
the waterfall phase and lead to copious production of monopole, monopole problem. One 
solution to evade the monopole problem is the shifted inflation model j|3] 



W = hS(-M 2 ^ nol-^M (3.45) 

String 



where and are defined in (|3.18|) and Mstring is a string scale. In this and next smooth 
cases, and have vevs during the inflation phase and monopoles are diluted away. However, 



in this section, we consider a simpler model, smooth hybrid inflation model ||5] in the context 
of the orbifold GUT model discussed in the previous section. For this purpose, we introduce 
a singlet chiral superfield S as before. Needless to say, this singlet field causes no change 
for the gauge coupling unification. Let us consider the superpotential for the smooth hybrid 
inflation 



Ws = S[-^ + ^-). (3.46) 



Here and are defined in ( |3.18| ) and we have omitted possible 0(1) coefficients. ( p.46| ) 
has U(1)r symmetry 

00 00 

S^e ia S, W->e ia W (3.47) 



under ( |2.201| ). SUSY vacuum conditions lead to non-zero VEVs for (0) = (0) = y/JIM, by 
which the PS symmetry is broken down to the SM one, and thus 



v PS = y/JM. (3.48) 



5 The renormalizable term, S(4>4>) 1 can be forbidden by introducing a discrete symmetry pq ], for example, 
— > —d> and (/>—></>. 
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U(1)r is conserved and no problem mentioned in subsection [2.9.71 occurs. Total Higgs su- 
perpotential including terms irrelevant to inflation is 



W = W PS + W s 



(3.49) 



where Wps is one given by (|3.21|) eliminating A15 term. | In the following analysis of inflation, 
we treat M as a free parameter with vp$ = \fJiM = 1.2 x 10 16 GeV fixed by the analysis 
of gauge coupling unification. Since M is involved in the non-renormalizable term, it is 
theoretically natural that M ~ M 5 ~ M GU t = 4.6 x 10 17 GeV Notwithstanding this 



condition is independent of the inflation scenario, we will find from the following analysis 
that M ~ Mqut is in fact consistent with the cosmological observations. Therefore, our 
GUT model is suitable for the inflation models with the parameters i>pg and Mqut fixed by 
the analysis for the gauge coupling unification. 

It would be worth mentioning that the Higgs superpotential includes a term W D H 6 ((f) 2 + 
2 ) through which all color triplets in and become heavy (( |3.21p ). Only the superpotential 
of Eq. (|3.46|) is relevant to inflation. 

Now we examine the smooth hybrid inflation scenario |35|]. The scalar potential from 
Eq. ( |3.46| ) is given byf\ 



V 



M 2 



+ 4S- 



M 4 



+ 



Considering the D-flatness condition, we normalize 



101 = 101 
and then V becomes 

V= (>- 
For a fixed a, V has a minimum at 



151 



x 1 



16M 2 



V2 



16M 4 ' 



X 



36a 4 + 16vt 



PS 



4^4 
3 a 2 ' 



(3.50) 



(3.51) 



(3.52) 



(3.53) 



where we have used an approximation for a 2 ^> v 2 PS (satisfied during inflation) in the last 
expression. The inflation trajectory is along this minimum and we obtain the potential along 
this path 



V ~ / 1 



27a 4 J 



(3.54) 



The potential diagram of the smooth hybrid inflation is depicted in figure 3.2 Note that 



6 If we accept more PS singlets other than S, we can construct an alternative PS invariant superpotential 
1. However it is less predictive for quark-lepton mass marices. 



7 Although supergravity effects from Kahler potential can play an important role |163 , we do not consider 
effects form supergravity in this paper, assuming a special Kahler potential for the inflaton field. 
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the PS symmetry is spontaneously broken anywhere on the inflation trajectory and thus, no 
topological defects such as strings, monopoles, or domain walls are produced at the end of 



inflation [H5 



The slow-roll parameters (e, if) and the parameter (£ 2 ), which enters the running of the 
spectral index, are defined as |36| 

r2 /T/'\ 2 QO A,f2„,8 A„A ^ 



M 2 P fV'\\ 32M>f s _ Av% s 



Tj — M 



V J 729(x 10 135(T 4 

., / V"\ 40Mj>v$s 
p 



V J 27a 6 



where the prime denotes derivative with respect to a. The slow-roll approximation is valid 
if the conditions, e and \r)\ <C 1, hold. In this case, the spectral index (n s ), the ratio of 
tensor-to-scalar fluctuations (r) and the running of the spectral index (a s ) are given by 

n s ~ 1 — 6e + 2r), 
r ~ 16e, 

a s = — -V ~ 16er] - 24e 2 - 2£ 2 . (3.56) 
dmk 

The number of e- folds after the comoving scale t§ = 27r//co has crosses the horizon is 
given by 

1 r ak V 9 

N k = — / da— ~ =— =- (a, 6 - a?) , (3.57) 

Mp J (Tf V 16M 2 4 S 1 fc /7 ' 1 ; 
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where ak is the value of the inflaton field when the scale corresponding to ko exits the 
horizon, and <7f is the value of the inflaton field when the inflation ends, which is determined 
by \rj\ = 1 so that 

40 



a) = ^M 2 P vU- (3.58) 
For a\ 3> cjj, several formulas given above are reduced into simpler forms, for example, 



f 27 

r 6 



N k - , n s ~ 1 - — — , a s ~ -77T72- ( 3 - 59 ) 



QVk 3AV 3iV, 



The number of e-folds iV& required for solving the horizon and flatness problems of the 
standard big bang cosmology, for k = 0.002 Mpc -1 , is given by [^6| 

where we have assumed a standard thermal history, and T r h is the reheating temperature after 
inflation. If gravitino has mass around 100 GeV as in the gravity mediated SUSY breaking, 
the reheating temperature is severely constrained (gravitino problem) from the fact that the 
gravitino decay products do not to destroy the light elements successfully synthesized during 
big bang nucleosynthesis ||164j| , 



T rh < 10 6 - 10 7 GeV. (3.61) 
The power spectrum of the primordial curvature perturbation at the scale k is given by 



1 V 3 / 2 27 



7V /2 — - ° Jl — (3 62) 



5 



This should satisfy the observed value by the WMAP f38j, V n ~ 2.457 x 10" 9 



Now we solve Eqs. (|3.57 ), ( |3.60|) and ( p.62|) and fix the model parameters of the inflation 



scenario. In our analysis, we have three independent free parameters, M, ak and T r h, with 
the fixed vps = 1.2 x 10 16 GeV. We solve the equations for a given T r h in the range 1 MeV 



< T r h < 10 GeV and find M and ak- Fig. [373] shows the results for M and ak as a function 



of T r h. We can check that a\ ^ a ^ and the slow-roll conditions are satisfied. The number of 
e-folds is depicted in Fig. [J7JJ Using these outputs and Eq. ( |3.55| ), we evaluate the spectral 
index (Fig. p.5|) , the tensor-to-scalar ratio and the running of the spectral index: 

0.963 < n s < 0.968, 
4.0 x 10~ 7 > r > 3.1 x 10~ 7 , 
-8.4 x 10" 4 < a s < -6.1 x 10" 4 (3.63) 

for 1 MeV < T r h < 10 7 GeV. The tensor-to-scalar ratio and the running of the spectral 
index are negligibly small. These results are consistent with the WMAP 5- year data p8| : 
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Figure 3.3: M and <Jk as a function of T r h. The solid and dashed lines correspond to M 
and <7fc, respectively. 




-2 2 4 6 8 10 



log 10 (T rh /GeV) 

Figure 3.4: The number of e-folds versus T^. 



n s = 0.9601^^, r < 0.24 (95% CL) and a s = -0.032+°;°^ (68% CL) (consistent with zero 
in 95% CL). The tiny value of r is in contrast with that of single field inflation, X a S a a = 
4,3,2,1,2/3, which gives rather large r around 0.1 (a = 2). In general, we do not 
need to impose either the left-right symmetry or M c = vps on the model. In this case, 
fps can be varied |53), and we repeat the same analysis for this general case with t>ps as a 



free parameter. Fixing, for example, T r h = 10 7 GeV, M can be obtained as a function of 
fps as shown in Figure |3.6| . Since M appears in the non-renormalizable term, it would be 
natural to identify M as an effective cutoff. Thus, the theoretical consistency leads to the 
condition, t> PS < M < M P , from which we obtain 1.9 x 10 14 GeV < v PS < 5.6 x 10 16 GeV. 



The number of e-folds and the tensor-to-scalar ratio versus vps are shown in Figs. [O] and 
|3~8| , respectively. 

The other outputs, the spectral index and its running, are found as 

0.967 < n s < 0.968, 
-6.4 x 10~ 4 < a s < -6.1 x 10~ 4 . (3.64) 

These are consistent with the WMAP data. 



70 



0.99 ; 

0.9 



0.96 ; 

0.95 ; 

0.94 ; 

0.93 : 

0.92 — 



2 4 6 8 10 

log 10 (T rh /GeV) 



Figure 3.5: The spectral index as a function of T^. The dashed lines correspond to the 
WMAP 5-year data, n s = OMO^H- 
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Figure 3.6: M as a function of fps (solid line). The lower and upper dashed lines specified 
the theoretical consist region for M, v-ps < M < Mp. 
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Figure 3.7: The number of e-folds versus vps- 
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Figure 3.8: The tensor-to-scalar ratio versus vps- 



3.5 SUSY Breaking Mechanism and Dark Matter 

For the scale fi > vps, there are only two independent gauge couplings a 4 and a 2 = a 2 R, 
so that the gauge coupling unification is easily realized with a suitable M c . Fig. |3l] shows 
the gauge coupling evolutions for \i > vps- For the scale \i > vps, there are only two 
independent gauge couplings a 4 and a 2 = «_r = «l, so that the gauge coupling unification 
is easily realized with a suitable M c . In this Fig. |3.9| , we have taken (corresponding to the 
result in the following neutral LSP arguments) 

M c = 2.47 x v PS = 2.95 x 10 16 GeV (3.65) 

and after including Kaluza-Klein threshold contributions into the gauge coupling evolutions, 
the gauge coupling unification is realized at 

M GUT = 7.54 x 10 16 GeV. (3.66) 

As M c is raised, Mqut becomes smaller. As mentioned before, we assume that a more 
fundamental SO(10) GUT theory takes place at Mqut- The origin of SUSY breaking and 
its mediation to the MSSM sector is still a prime question in any phenomenological SUSY 
models. Since flavor- dependent soft SUSY breaking masses are severely constrained by the 
current experiments, a mechanism which naturally transmits SUSY breaking in a flavor-blind 
way is the most favorable one. 



In higher dimensional models, the sequestering |165| is the easiest way to suppress flavor- 
dependent SUSY breaking effects to the MSSM matter sector. Since all matters reside on the 
PS brane in our model, the sequestering scenario is automatically realized when we simply 
assume a SUSY breaking sector on the brane at y = 0. SUSY breaking is propagated to PS 
brane by the SM gauge multiplet in the bulk (gaugino mediation) ||130 |. There is another 



SUSY breaking scenario. That is, gauge multiplets are also confined on MSSM brane and 
SUSY breaking is transmitted only by supergravity effect, the soft terms on the visible brane 
can be understood by the anomalous violation of a local superconformal invariance and called 
anomalous SUSY breaking (AMSB) ||165|| , [ 132j| . AMSB suffers from the problem that the 



slepton masses become tachyonic. There are positive contributions from modulus [|166|| and 
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Figure 3.9: Gauge coupling unification for M c > vps- Each line from top to bottom 
corresponds to g 3 , g 2 and g\ for fi < vps, while g 3 = g± and g 2 = g 2 R for fi > vps- Here, we 
have taken M c = 2.47 x vps- 



D-term. It should be remarked that SO(10) GUT including high dimensional representation 
discussed here have U(1)b-l an d D term. So in general we can combine modulus and D- 
term contributions ||167| . Anyhow we adopt the gauge multiplet in bulk and do not discuss 
AMSB. 

Gaugino mediation The SO(10) gauge multiplet is in the bulk and can directly com- 
municate with the SUSY breaking sector through the higher dimensional operator of the 
form, 

£ ~ 8{y) ! d 2 9 ^tr [W a W a ] . (3.67) 

•J 5 

Here X is a singlet chiral superfield which breaks SUSY by its F-component VEV, X = 8 2 Fx- 
Also W is the field strength chiral superfield constructed from the vector superfield V and 
chiral covariant D a , 

W a (x,6,6) = --D+Dt (e- v D a e v ) 

= -i\ a {x + lOa^O) + [5£D{x + iQo$) - i^^F^x + iOa^)] 9 P 

+ Oeia^dPXjiix + iOajS) (3.68) 

and 

D a = ^~ ^{^) a ^,. (3.69) 
Therefore, the bulk gaugino obtains the SUSY breaking soft mass, 

F X M C F x ( 



M\ ~ - — [ — (3.70) 

M| M P \MpJ 
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where M c comes from the wave function normalization of the bulk gaugino, and we have 
used the relation between the 4D and 5D Planck scales, M|/M c ~ Mp in the last equal- 
ity. As usual, we take M\ =100 GeV-1 TeV. Once the gaugino obtains non-zero mass, 
SUSY breaking terms for sfermions are automatically generated through the RGE from the 
compactification scale to the electroweak scale. This scenario is nothing but the gaugino 
mediation |pL 30|| and flavor-blind sfermion masses are generated through the gauge interac- 



tions. In this setup, a typical gaugino mass in Eq. ( 3.70 ) is smaller than the gravitino mass 
m 3/2 ^ F x /Mp by a factor M 5 /M P < 1. 

Vps, we find that the right-handed slepton (normally, 



As discussed in Ref. 34 



for M c 

stau) is the LSP, because the sfermion mass spectrum is obtained from the boundary condi- 
tion with vanishing soft masses at M c = vps = 1.19 x 10 16 GeV. This result is problematic 
for cosmology. As pointed out in Ref. [ 168 ], this stau LSP problem is cured by the soft mass 
RGE running from the compactification scale to the GUT scale in a GUT model. In the 
following, we apply this idea to our model with M c > vps- 

For the scale, vps < t 1 < M c , we are in the PS stage and the RGEs of gaugino and 
sfermion masses are given by 



d_ f Mil 
dt V a 2L 



d /M 4 
dt \ «4 

dm\ 15 _ 3 „, 2 

dt 47r 2tt 



d_ (M 2R 
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&2R 
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47T 



a 4 M 4 



2tt 



a 2R M' : 



2_R) 



(3.71) 



where t = ln(/z/M c ), a 4 and ck 2 l = «2_r are the PS gauge coupling of the corresponding 
gauge groups (whose RGE solutions are obtained in the previous section), and M 4 , M 2L and 
M 2 r are the corresponding gaugino masses. Sfermion mass spectrum is obtained by solving 
the RGEs with the boundary conditions, m F (M c ) = m Fc (M c ) = 0. Analytic solutions of 
Eq. (|3.71|) at fi = vps are easily found: 
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«2L(M C ) V 1 
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<*2fl(M c ) Y 
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.(3.72) 



Note that the PS model is unified into a more fundamental SO(10) model and this unification 
leads to the well-known relation, 



M 4 

Q! 4 



M- 



2L 



M- 



2R 



M l/2 



(3.73) 



«2L &2R «GUT 

where M\/ 2 is the universal gaugino mass at the unification scale. Thus, the formula as for 
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sfermion masses are simplified as 



m 2 Jv PS ) = m 2 p, c (vp S ) 



/M 1/2 V 

V«GUT/ 



(a 2 4 (M c 



- <4(v PS )) + - «(M C ) - a 2 2L {v PS )) 



(3.74) 



Solving the RGEs in the MSSM with the universal boundary condition m F (vps) = 
m 2 p c {vps) — m o at /i = vps, we obtain the sfermion masses at the electroweak scale. In 
our model, the PS scale is almost the same as the usual SUSY GUT scale in the MSSM 
(Mqut — 2 x 10 16 GeV) and the gauge couplings are roughly unified at the PS scale, 
on{vps) — a 2 (v PS ) = a 2 p (see Fig. gjg) . 




Figure 3.10: The ratio, R = m 2 c /M 2 26 , as a function of M c /vps- Here, soft masses have 
been evaluated at \i = Msusy- R = 1 when M c /vps = 2.47. 



Therefore, our study on the sfermion masses are almost the same as the one usual in 
the constrained MSSM. For a small tan/3 (say, tan/3 = 10), we neglect Yukawa coupling 
contributions to the soft masses of right-handed sleptons, and the analytic solutions of the 
MSSM RGEs are found to be 

M 126 (/i) = a^fMlll), 

V«GUT/ 

™'M = f— ) £r[<*l(vps) ~<* 2 M] +ml (3.75) 
\ o;gut / 1 1 

If m is large enough, the slepton (stau) mass is bigger than the bino mass (Mi 26 ). In our 
model, mo is given as a function of M c . Fig. p.lOj shows the ratio, 



2 

m 



R^TT^i (3-76) 



Mf 26 



as a function of M c /vps- We can obtain R> 1 for M c /vps > 2.47. 

Now, for M c /vps > 2.47, the bino-like neutralino will be the LSP and a good candidate 



for the cold dark matter in cosmology |169| . For a small tan/3, the annihilation processes of 
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the bino-like neutralino are dominated by p-wave and are not so efficient. As a result, the 
neutralino relic density tends to exceed the upper bound of the observed dark matter density. 
This problem can be avoided if the neutralino is quasi-degenerate with the next LSP slepton 
(stau), and the co-annihilation process with the next LSP can lead to the right dark matter 
relic density. In our result, such a situation appears for M c ~ 2.47 x vps — 2.95 x 10 16 GeV. 
It would be interesting to note that the discrepancy of the abundance of 7 Li between the 
observed values in WMAP and in metal poor halo stars may be explained the degeneracy 
between the LSP neutralino and stau [170]. 



3.6 Leptogenesis in 5D 



In section |2.6j , we discussed non-therml leptogenesis in 4D, where A = 1(T 8 . This gave rise 
to a naturaliness problem. In the following, we discuss how this problem is modified in the 
smooth hybrid inflation model |. 

It goes from ( |3.19| ) and ( |3.46| ) that the interaction terms relevant to inflaton decay are 



The potential V is 



W = S 



V 



-/i 2 + 



s 



M 2 



+ 



Y v 

R rpc fpc 



+ 



R 77>c fpc 



M 2 ' Mr^ Ri Rj 



(3.77) 



(3.78) 



So SF^Fftj interaction term is 



M PS ~c ~c 

Hspp = 8 M 2 M 5 YRSFRiFRj '' 



and the decay ratio at tree level is 
Thus A corresponds in this paper to 



A 2 « Y? 



64 2 M* s 1 
8tt r M 2 M a M s ' 



R MqM s ' 



(3.79) 



(3.80) 



(3.81) 



8M|, 



where we assumed M = M 5 . M s is obtained from ( |3.78| ), M| = 01 ^[i s ■ Substituting the 

values obtained, we get Y R = 1.0 x 10" 6 if BR{S -> F c Rl F c R1 ) = 1. 
In the SO(10) in 4D model, 



M m = 1.64 x 10 11 , Mr2 = 2.50 x 10^, M R3 = 8.22 x 10 lz GeV 



12 



,12/ 



(3.82) 



8 In ths section we discuss on the sneutrino case (S ee se ction |2.6| for neutrino.). Also there is a contribution 
from soft SUSY breaking terms (soft Leptogenesis) [171|. 
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In this model Ymj was related with quark- lepton masses via 126 Higgs field. On the other 
hand, as is well known from (|3.19| ) Y Pi j in the present model are independent on the quark- 
lepton mass matrices. However, Dirac neutrino mass matrices Y ui j are same in both models 
and therfore eigen values of Mr may not be so different to each other. So M s > M R3 , ikL 



and all decay channels mentioned above are open with the same decay ratio from SUSY |84 

Single inflaton model suffered from the fine tuning problem A = O(10~ 14 ), which was one 
of the motivation to have introduced hybrid model. As we mentioned above, this problem 
has reappeared as Yr = 1.0 x 1CT 6 . This constraint comes from Tr < 10 6 GeV. However, 
we are free from this condition if gravitino mass is of order 100 TeV. Let us consider ( [3.67 ) 
in more details, and we put 

C = c g 6(y) I d 2 6 -^tr [W a W a ] . (3.83) 
From which we obtain the gaugino mass 

= c ^{w P ) (3 - 84) 

If we adopt c g = 0.1, then we obtain ~ lOOTeV heavy enough to be free from gravitino 
problem and from the constraint Tr < 10 6 GeV. One gravitino decays to one LSP and we 
obtain 

n L sp/s = n 3/2 /s. (3.85) 
LSP as a DM must satisfy the observational constraint, 

^Lsph 2 < 0.1, (3.86) 

which leads us to 

n 3/2 /s = n L sp/s = 3 x 10~ 12 x (mLSp/lOOGeV)- 1 x (n LSP h 2 /0.1). (3.87) 

If gravitino is generated in the reheating heat bath after the inflation, 713/2 and Tr are related 
with 

n 3/2 /s ~ 1.5 x 10~ 12 x (T R /W 10 GeV) (3.88) 

and therefore 

T R <2x 10 10 GeV x (mLSp/WOGeV)- 1 x (n LSP h 2 /0.1) (3.89) 
So Mr must be non-thermally generated. 
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3.7 LFV in 5D 



Recently MEG collaboration [|1 72|| reported new results of a search for the /i — > cy decay 
and a maximally likelihood analysis sets an upper limit at 90% C.L. on the branching ratio, 

BR(/i ->• erj) < 2.4 x 10" 12 , (3.90) 

Before the data of (|3.90| ), we heard that they found 3 events as the best value for the 



number of signals in the maximally likelihood fit, which corresponds to [|173 

BR(/i ->■ erf) = 3 x 10~ 12 



(3.91) 



for the center value. This result was denied and reduced to ( [3.901 ) • However, our model gives 
rather large branching ratio marginal to the experimental upper bound (see Fig. |2.3| ). So it 
is very interesting to consider what happens if the observation gives non-null result around 



this value 11174 



The relic abundance of the cold dark matter (CDM) in 2a range has been measured as 

n CDM h 2 = 0.1131 + ±0.0034. (3.92) 

As is well-known, in SUSY models with the i?-parity conservation, a neutralino, if it is the 
lightest sparticle (LSP), is the promising candidate for the CDM in the present universe. 

For tan (5 = 45, \i > and Aq = 0, for example, we can find more precise relation between 
m Q and My 2 than (|2.89|) such as 



m (GeV) = 125.3 + 0.329M 1/2 (GeV) + 5 x 10~ 5 (M^GeV)) 1 



(3.93) 



along which the observed relic abundance, flcDMh 2 = 0.113, is realized. In our analysis, 



we have employed the SoftSUSY3. 1 .4 package ||175|| to solve the MSSM RGEs and produce 
mass spectrum. Then, the relic abundance of the neutralino dark matter is calculated by 
using the micrOMEGAs 2.4 |176| with the output of SoftSUSY in the SLHA format [|i~77f . 

Anyhow let us continue our arguments without adopting ( |3.91| ). Through the seesaw 
mechanism M, the light neutrino mass matrix is given by 



m v = Y u T M^Y u v 2 u 



UmnsDMmns 



(3.94) 



in the basis where the mass matrix of charged lepton is diagonal with real and positive 
eigenvalues. Here v u is the VEV of the up-type Higgs doublet in the MSSM, D u is the 
diagonal mass matrix of light neutrinos, and Umns is neutrino mixing matrix. This is 
equivalent to the expression of the right-handed neutrino mass matrix as 



M, 



Y U U TBM D U 1 U^ BM Y l/ 



T 



(3.95) 



Once the information of the Dirac Yukawa coupling, the mass spectrum of the light neutrinos, 
and the neutrino mixing matrix is obtained, we can fix the right-handed neutrino mass 



matrix. In order to determine Mr, we follow the manner in [47 . 
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We consider the normal hierarchical case for the light neutrino mass spectrum, for sim- 
plicity, and describe D v in terms of the lightest mass eigenvalue m,\ and the mass squared 
differences: 



D v = diag mi 



Amf 2 + m\ , 



Amj 3 



Here we adopted the neutrino oscillation data [|178| : 



Am 2 12 



7.59 x 10" 5 eV J , Am( 3 = 2.43 x lO -3 eV 2 



(3.96) 



(3.97) 



Contrary to the case of the minimal SO(10) GUT in 4D, right-handed heavy neutrino mass 
matrix is independent on the charged fermion mass matrices, and we assume the neutrino 



mixing matrix of the so-called tri-bimaximal form [|179 



U- 



TBM 



( 



\ 



(3.98) 



which is in very good agreement with the current best fit values of the neutrino oscillation 
data [178]. | As has been discussed above, the data fit for the realistic charged fermion mass 



matrices is same as in the minimal SO(10) model, and we here use the numerical value of Y v 
in Eq. ( [2.85 ) at the GUT scale for tan/3 = 45, in the basis where the charged lepton mass 
matrix is diagonal. Substituting Eqs. ( gjgg ), ( gg§ ), (|337Q and ( gggp into Eq. ( ggg ), we 
obtain the right-handed neutrino mass matrix as a function of only mi. Since it has been 
shown [f|7]] that the simple 5D SO(10) model can reproduce the observed baryon asymmetry 
of the present universe for mi ~ 1.8 x 1CT 3 eV through non-thermal leptogenesis, we take 
this for a reference value of mi. In this way, M R is now completely determined, but not 
yet diagonalized. Changing the basis where M R is also diagonal, we find the neutrino Dirac 
Yukawa matrix, 



-0.00119 + 0.00002682 
0.00135 + 0.001672 
-0.0265 -0.01732 



-0.00108 - 0.0004852 
-0.0253 + 0.00154? 
0.0609 + 0.0275? 



-0.000392 + 0.000421i 
0.0237 + 0.0008512 
0.790 - 0.04362' 



,(3.99) 



and the diagonal right-handed neutrino mass matrix with eigenvalues (in GeV) 

M Rl = 1.03 x 10 10 , M R2 = 7.55 x 10 11 , M Ra = 3.22 x 10 15 . (3.100) 



l R 2 



l R 3 



Reflecting the independence of M R from the charged leptons, Y v of ( |3.99| ) is different 
from ( |2.85 ). Now we are ready to analyze the // — > decay rate by using the completely 
determined neutrino Dirac Yukawa coupling matrix and the right-handed neutrino mass 



This is found now correct upto a zero-th order approximation after the non zero discovery of #13 (1.3). 
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7.0 x 10-' 2 




M1/2 (GeV) 



Figure 3.11: The branching ratio BR(/i — > e'j) for the minimal SO(10) model as a func- 
tion of Mi/ 2 (GeV) along the cosmological constraint of Eq. ( |3.93|) . The short dashed line 
corresponds to the MEG result, BR(/i — > ej) < 2.4 x 1CT 12 . This diagram is same as the 
left-upper panel of Fig. 2. 3 except for the use of Eq. ( |2.89| for the latter case. 



eigenvalues^, Eqs. (|2.84j) and ( p.85| ) for the minimal SO(10) model while Eqs. (|3.99|) and 



( p.l00| ) for the simple 5D SO(10) model. For tan/3 = 45, A = and fi > 0, the branching 
ratio of — > for the minimal SO(10) model is shown in Fig. 3.11 as a function of the 
universal gaugino mass at the GUT scale, along the relation of Eq. ( p. 93 ) to satisfy the 



observed relic density for the neutralino dark matter. The short dashed line corresponds to 
the branching ratio of Eq. ( |3.90| ). The universal gaugino mass becomes Mi/ 2 > 840 GeV. 



Fig. 3.12 depicts the same result as in Fig. 3.11, but for the simple 5D SO(10) model. 
In this case, the universal gaugino mass is relatively low and becomes M1/2 > 440 GeV. 
This is because the components in the neutrino Dirac Yukawa matrix relevant to the LFV 
between 1st and 2nd generations are smaller than those in the minimal SO(10) model and 
a lighter sparticle mass spectrum is necessary in order to achieve the same branching ratio 
(see Eq. flTTp). 

The set of the CMSSM parameters proposes a good benchmark point for the SUSY search 
at the LHC and we presented (sparticle) mass spectra for each model in Table 3.3, along 
with other observables which can be compared with the experimental bounds [ 174 | before 



the dicovery of the Higgs-like particle. (We remained the old data before the anouncement at 
July 2012 by the LHC since they are useful for exhibiting what should be correcyted. We will 
discuss on the impact of the recent LHC results at the end of this section.) In both results, 
the lower bound on the Higgs boson mass rrih > 114.4 GeV [|180|1 is satisfied. However, 
squark masses in 5D SO(10) model are disfavored. Other phenomenological constraints we 



in 



Previous analysis for the minimal SO(10) model with general parameter sets, see [BOj 



SO 



Table 3.3: Mass spectra (in GeV) and phenomenological constraints for the two SUSY 
SO (10) models with the universal boundary conditions in the CMSSM cited from ||174|| . A s 
in the first and second frames are the trilinear term and CP-odd neutral scalar, repectively. 



See also Table 1373 reflecting the discovery of Higgs-like boson around 125 GeV by the LHC. 
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Figure 3.12: The same figure as Fig. 3.11 but for the simple 5D SO(10) model. 



consider here are 



2.85 x 1(T 4 < BR(6 -> 57) < 4.24 x 10" 4 (2a) 
BR(B S n+fT) < 5.8 x 1(T 8 



181 



1821. 



(3.101) 
(3.102) 



Q The results for the minimal SO(10) model satisfy these constraints, while BR(6 — > 57) for 
the simple 5D SO(10) model is marginal (about 3.4a away from the center value). 

We can accommodate our results with muon anomalous magnetic moment, another im- 
portant signal of New Physics BSM. 

Table 3.3 also includes the sparticle contributions to Aa^. We can see that the result 
for the minimal 5D SO(10) model favors the deviation obtained by using the data of the 
hadronic r decay, while the deviation from the e + e~ data is favored by the result for the 
simple SO(10) model. These are the results before the anouncement of the LHC at 7 Tev 
and 8 TeV run. 



GUT and the LHC results 

We study the impact of the recent LHC discoveries of Higgs-like object around = 125 
GeV, h — > 77, SUSY search etc. in connection with SO(10) GUT models. Here the consis- 
tencies of the SUSY have been analyzed without considering the details of model but with 
the pattern of soft SUSY breaking. 
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( P.102| ) is updated to BR(B S < 4.5(3.8) x 1CT 9 at 95% (90%) C.L. [|L83 
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One loop correction to the lightest Higgs mass in CMSSM is ||184| [|185 



3 TTt^ 

m 2 ^M 2 cos 2 2/3+— 



In — % 



+ 



mi 



Ml 



1 - 



YIMl 



(3.103) 



where 

M s = y/m h m h , X t = A t - //cot/3, v = 174GeV (3.104) 

with the trilinear Higgs-stop coupling constant A t . For large tan/3, some negative corrections 
appear 

h\v 2 /i 4 



Am 2 h 



from sbottom, and 



Ami 



16vr 2 Ml 



h\v 2 fi 4 



48vr 2 m% 

from stau. Here the bottom Yukawa and the tau Yukawa are 



h h 



m b 



fcos/3(l + tan/3A/ifc 



h r 



m T 



vcos/3(l + tan/3A/i T ) 



(3.105) 
(3.106) 

(3.107) 



with one-loop corrections of Ah^ and Ah T . Anyhow these corrections may be subdominant. 
The recent review of SUSY search at the LHC gives very severe constraint on CMSSM 
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Here we set Higgs-like object around 125 GeV as lightest Higgs h, and (|3.103|) indicates 
(a) rather large stop mass or (b) large A t . (a) implies a large discrepancy between m t and 
mi and is in the inverse direction of that of original SUSY motivation, the loop cancellation 
in Higgs mass hierarchy. This large stop mass need large m or large M1/2 for gaugino 
mediation, which is very severe from the LHC [ 186|| . We did not adopt the second choice 
since we set Aq = 0. We will be back on this point later. 

As a result, someone assert that CMSSM is strogly disfavored ||187|| . If it is indeed 



the case, it should be considered very worrisome since our analyses of LFV, leptogenesis, 
sparticles mass spectra have been based on the universal boundary condition of CMSSM 
from (|2.80| ) to fl2.83| ). As we mentioned, these universal boundary conditions are natural 
except for ( |2.81| ) in the framework of SO(10) group. It is more natural that we relax ( [2.81| ) 
and take mu u and mn d as free parameters (NUHM1 or NUHM2). It is not obvious to us 
that it improves the problems [ISS]. Here we only point out that this may suppress sfermion 
masses since 

SC a (t)g 2 a \M a \ 2 + ^Y ig 2 S, (3.108 

a=l,2,3 ° 



dt ^ 



where 



S = Tr[Yjm z 



m 



mi + Tr[m 2 - m 2 - 2m| + m§ + m 2 



(3.109) 



Gaugino mediation which we adopted in 5D model makes the situation worse. This is because 
all matters get their masses from gaugino by RGE and gaugino mass must be unacceptably 
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large. Also we may generalize Q2.82Q : if the hidden sector field is not the singlet unlike ( |3.67D 
but, for instance, 54-plet 

£ = J d 2 6j^W a W? } (3.110) 

we have non- universal gaugino masses at GUT scalej lSS ] 

M 3 : M 2 : M x = 2 : -3 : -1 (3.111) 



in place of ( 2.82 ). So, even if we adopt CMSSM or minimally extend it, we must first 
construct the most suitable soft mass conditions (m , Ml/2, A , tan/3 for CMSSM or 
mo, Mi/2, A , tan/3, /i, for NUHM2 ) without singlet and with singlet (the next- 
to-minimal supersymmetric model (NMSSM) ||190|| and the nearly minimal supersymmetric 
model (nMSSM) ||191|| etc.) due to several SUSY breaking mechanism. 

The NMSSM and nMSSM are primarily motivated by // problem and introduce gauge 
singlet superfield S. 

The superpotential of NMSSM is 

W = Wmssm + XSH u H d + ^kS 3 , (3.112) 

where Wmssm is the superpotential of the MSSM without /i term. The F term is 

V F = \X\ 2 \S\ 2 (hIH u + H\H^ + \kS 2 \ 2 (3.113) 

S couples with Higgs only through W and has no gauge coupling. Therefore, D-term is same 
as MSSM 

v D = \gl\HlH d \ 2 + \{gl + gl) {h ] u h u - H\H d ) 2 . (3.114) 

V S oft = m 2 H E\E u + m 2 H H\H d + m 2 s \S\ 2 + (\A x H u H d S + ^A K S 3 + c.c)j . (3.115) 
The mass of the lightest Higgs boson becomes at tree level 

K MSSM ) 2 < m 2 z (cos 2 m + 2|A|2 2 S f , (3.116) 

V <?L "T ^2 / 

though it works well for small tan/3 around 1. 

The nMSSM is a model whose superpotential is 

W = W M ssm + ZfM 2 S. (3.117) 

In this model DM is light singlino dominated neutralino, which greatly enhance the decay 
of h into light neutralino. This greatly suppresses Br(h — > 77) in comparion with the SM 
and is incompatible with the LHC observation. 

One of the peculiar properties of our data fiiting among others are large tan/3. 



84 



It goes from (|2.87 ) that this A t ~ A T m Aq value enhances LFV in general. 



ml 



m| 3 + A-gL v (A* T cos(3 - ny T smj3) 

v(A T cos(3 — fi*y T sinf3) rrie 3 + Ag fl 



where A is the contribution from D term and 

A^ = (T30 - Q 4> sm 2 9 w )cs(2f3)m 2 z . 
For large tan/3 value, the following value is a check point [|192| [|193 

6 



(3.118) 



(3.119) 



x 



3.5 x 10 



m. 



-5 



tan/3 
50 



l.hps 



230MeV 



\V f 



e//| 



0.040 



[l6n 2 ) 2 e 2 Y 



Here 



with 



x 



(1 + e 3 tan/3) 2 (l + eotan/3) 1 



1 A-^cotp Q/ uf etc _ 



16tt 2 



— , H 2 (x,y) 



xlnx 



+ 



y\ny 



(3.120) 



(3.121) 



(3.122) 



H' ■ [l-x)(x-y) (l-y)(y-x)' 

There may be 3a deviation from the SM in h — > 77 [|194j| . Its decay ratio is given by ||195 l 



+ ^ ± <4^ ± ) + ^4 ± 44^ ± ) + ^iV c/ e^(r / ) I (3.123) 



where r a = imy L j2m a ) 2 . Ai (i specifies spin here) are defined by 

A x {r) = -[2r 2 + 3r + 3(2T-l)g(r)]/T 2 
A 1/2 (t) = 2[r+(r-l)g(r)]/r 2 
Mr) = -[r-g(r)]/r 2 , 



where 



arcsm Wr 



9{r) 



log v , — 27T 



for t < 1 



for r > 1 



(3.124) 
(3.125) 
(3.126) 



(3.127) 



is negative definite, and Aq and A1/2 are monotonicly increasing positive definite func- 
tions. Thus W boson loop dominates and top quark suppreeses the decay rate in the SM. 
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However, if the contributions of sfermions dominate over the W boson loop, then T(h — > 77) 
may be enhanced in comparison with the SM. Taking all these situations into consideration, 
we will consider how these data constrain the model. These properties are fit to the LHC 



data if we select parameter regions, for instance, as ||185 | 

A T (trilinear) « 500 [GeV], fj, « ITeV tan/3 « 60, (3.128) 

which gives 

Br(h -»■ 77) « 1.5Br(h -> 77) 5M . (3.129) 

In this case a lighter stau mass m? = 135 GeV is very near to rrih- They give another 
example 

A T m 1500 [GeV], fi « 1030GeV tan/3 « 60, (3.130) 

which gives a lighter stau mass m f = 106GeV < m^. 

Using Sof tSUSY3 . 1 . 6, we give an example in the table [3]4| for tan/3 = 45 realized in 
the minimal SO(10) GUT in 4D 0. More elaborate sample calculations have been made 



in the CMSSM and NUHM1 |88| and in the pMSSM p7fl . However, SUSY GUT has 
its advantages over MSSM and pMSSM etc. since it can fix the Yukawa coupling and 
tan/3. LFV and leptogenesis processes need the detailed structure of the Dirac Yukawa 



coupling Y v . As you can see from Figure pT3| , M l / 2 must satisfy 350 < Mi/ 2 < 800 [GeV] if 
10 < <5a M x 10 10 < 40 in the case of trilinear A Q = 0. This is beacause 5a ^ monotonically 
decreases as Mxii increases. 

However, it can be relaxed for A ^ and larger M1/2 may be accepted. 

After almost completing this review, LHCb Collaboration has anounced that they have 
found the first evidence of Br{B s -> = (3.2±Jj) x 10~ 9 pg| , which is in agreement 

with the SM prediction, Br(B s = 3.23 ± 0.27) x 10~ 9 . 



Discussion 

We have reviewed SO(10) GUT mainly based on our works. The most peculiar property 
of minimal SO(10) GUT is its high predictivity and its predictions cover over all particle 
phenomena and over a wide range of energies, from 10 -5 eV to 10 16 GeV. This is a rather pe- 
culiar property of the minimal SO(10) model. This high predictivity, in other words, implies 
that this model suffers from continuous hard checks from many observations. Furthermore, 
not only data fitting over a wide range but also the internal consistencies like the compati- 
bilities with a gauge coupling unifications and with the detailed SUSY breaking mechanism 
etc. have been considered on the minimal SO(10) GUT. 

We first discussed the minimal SO(10) GUT in 4D. The data fitting with quark masses 
(6), CKM mixing angles and CP phase (3+1), charged lepton masses (3) were fitted with 
14 parameters (see section 2.4.2 for a more generic form). Adjusting remaining only one 
parameter and one scaling factor cr, we could fix full neutrino oscillation data; mass square 

12 We are grateful to H.Ishida for the help of numerical calculations. 
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Table 3.4: Mass spectra (in GeV) in various universal boundary conditions which give 
around 125 GeV. In all columns except for the first one we used (|3.93|) . 





mass 


in units of GeV 




' ' U 


2000 




1000 


1300 




1600 




M1/2 


2000 




2000 


2500 




3000 




A 


-3500 


-1800 


-2700 


-3500 


tan/3 


45 




45 


45 




45 




m ho 


125.1 




123.7 


124.7 




125.5 




Ho 


2266 




2011 


2504 




2982 




A 


2266 




2011 


2504 




2982 




H± 


2267 




2013 


2505 




2984 




Trig 


4268 




4218 


5204 




6179 






8702 


1656 


8615 , 1639 


1085. 


2059 


1309. 


2478 




-2754 , 2757 


-2422 , 2426 


-3057 , 3060 


-3663,3666 




1656 


2757 


1640 , 2426 


2059. 


3060 


2478. 


3666 




2757 


2422 


3058 


3665 


Ul,2 


4280 


4134 


3926 , 3762 


4849. 


4643 


5762. 


5513 


Cl,2 


4280 


4134 


3925 , 3762 


4849, 


4643 


5761. 


5513 


tl,2 


2895 


3458 


2917,3368 


3548, 


4113 


4179. 


4856 


di, 2 


4281 


4118 


3926 , 3743 


4849, 


4619 


5762. 


5484 


h,2 


4281 


4117 


3926 , 3743 


4849, 


4619 


5762. 


5483 


h,2 


3428 


3560 


3315,3383 


4064. 


4134 


4811. 


4881 


ei, 2 


2390 


2133 


1655 , 1248 


2095. 


1599 


2534. 


1951 


#1,2 


2389 


2130 


1655 , 1247 


2094, 


1597 


2533, 


1948 


n,2 


1132 


2027 


6061 , 1476 


6609. 


1834 


7410. 


2197 


u el,2 


2388 




1653 


2093 




2532 






2387 




1653 


2092 




2531 




V T 1,2 


2021 




1466 


1825 




2190 




Bt(B s -)• fi+fj,-) 


3.9 x 10~ 9 


7.6 x 10~ 9 


2.9 x lO" 9 


6.2 x 10~ 10 
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differences (2), MNS mixing matrix (angles (3), CP phases (3)) as well as the Dirac neutrino 
(Y v ) and heavy right-handed neutrino (Mr) mass matrices. Y v and Mr are very essential for 
LFV and Leptogenesis, respectively. Unfortunately our predictions were found to deviate 
after Kamland in # 13 and the mass ratio Am^/Am^. However, this does not mean that 
minimal SO(10) is excluded. Our parameter search is natural but not exhaustive for the 
following reasons. 

1. We set a, (3 (see section 2.4.2) fixed to or it for simplicity. 

2. We assumed only type I seesaw. Minimal SO (10) model allows both type I and Type 

II. 



The general case was discussed in pq| . Their data fitting was improved but not in very 
good agreement yet. More important is that it suffers from some internal problems like the 
instability of gauge coupling unification due to intermediate energy scales and the no-go 
theorem on the SUSY breaking etc. 

3. If we relax the "minimal" condition but preserving renormalizability, we may add 120- 
plet Higgs. This may improve the conflict between the gauge couplin unification and several 
intermediate energy scales. Unfortunately, no one has succeeded in it explicitly. 

Both on the detailed mismaches of data fitting and on the internal consistency like the no- 
go theorem of U(1)r symmetry, there are several loopholes, making us stay in 4D. However, 
it seems a natural choice for us to go to extra dimensions. This is rather reasonable since 
our bottom-up approach should have a harmonic encounter with the top-down approach at 
GUT scales. There it is very natural for GUT to have some remnants of extra dimensions. 
Thus we have proceeded to 5D SO(10) orbifold GUT. 

Dimension five may not be our final goal but it must provide essential merits for extra 
dimensions if it is real. 

We have repeated a comprehensive analyses in 5D which was done in 4D and showed how 
they improve the situations compared with those in 4D. The problems and deficits discussed 
in section 2.8 are remedied in 5D. 

4. Data fitting is expected to be improved since Mr is free from M126 unlike the original 
SO(10) case. 

5. Superpotential (|3.49|) breaks the PS symmetry to the SM group without any interme- 
diate energy scale and does U{1)r consistently. 

Though we have accepted unrenormalizable terms but they are strictly guided by the 
renormalizable action and inherit the merits of the renormalizable model. 

The recent LHC discoveries of a Higgs-like boson around 125 GeV and of the null result 
of SUSY particle searches etc. give very important constraints on the SO(10) GUT model 
in 4D and 5D, which forces us to make some modifications to them. Neverthless we have 
cited the old results some of which are not compatible with the LHC results. This is because 
it is very important to first clearify the discrepancies between the theoritical predictions 
and observations. On going and near future experiments of accelerators and cosmological 
observations must drive GUT to precision science. 
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Appendix A 

SO(10) Group Properties 



A.l Spinorial representation 

This section is based on [^TJ 

SO(2N) is expanded in SU(N) basis, Xi (* = 1> ■■■■>N), constitutes Clifford algebra, 

{Xii Xj\ &ij 

{K,Xj} = {xlx]} = (A.l) 
Now let us define the 2N operators which constitutes vector basis, 

r 2j-i = -i{Xj - x)) 

T 2J = (Xj + x}) j = l,2,...,N (A.2) 

satisfies 

{r M , T u } = 25^ u fx, u = l,...,2N (A.3) 
Then we can construct a representation of SO(2N), where 

E^ = i[r^, r„] (A.4) 

are the 45 dimensional generators. The spinor representation in SO(2N) is 2 N dimensional 
and is split into 2 N ~ 1 dimensional representation under a chiral projection operator, 

r = z JV r 1 r 2 ...r JW (A.5) 

For the case of SO(10), 

|16_> = 116*) = xlxlxlxlxtm + ^ klm xlxlxlm% + xH0>S* 
|16+> = \16 R ) = |0)1 + \xlx]\0)^ j + Y/ 3klm x)xlx\xl\m^ (A.6) 
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The Yukawa coupling of chiral Fermions with Higgs scalars is given by 

#rv..r p U'$ / ,.. p (A.7) 

Here <& M ... P is Higgs scalar in tensorial representation which is discussed in the next section. 

A. 2 Tensorial representation 



The arguments of this appendix is based on |JD| . Our formulation is in tensorial representa- 
tion. 

Here we first introduce Y diagonal basis: 1 = V + 2'i, 2 = V -2'i, 3 = 3'+4'i, 4 = 3'-4'i, 
5 = 5' + 6'i, 6 = 5'- &i, 7 = 7' + 8% 8 = 7'- 8% 9 = 9' + 0% = 9'- 0% up to the 
normalization factor It is more convenient since (1,3,5,7,9) transforms as 5-plet and 

(2, 4, 6, 8, 0) transforms as 5-plet under SU(5) x U (l)x (for that reason Y diagonal basis could 
also be called SU(5) basis). Consequently, (1, 3) and (2, 4) are SU(2) L doublets with definite 
hypercharges Y = | and Y = — |, respectively. Similarly, (5,7,9) and (6,8,0) transform 
under SU(3)c as 3 and 3 with definite hypercharges Y = — | and Y = |, respectively. 

Note that under the complex conjugation (c.c), 1 = 2, 3 = 4, 5 = 6, 7 = 8, 9 = 0, and 
vice versa. The 50(10) invariants are build in such a way that an index a is contracted 
(summed) with the corresponding c.c. index a, for example, T... a ...T.. ...... 

The basis in A = 45, £> = 120, $ = 210 and A + A = 126 + 126 dimensional spaces 
are defined by totally antisymmetric (unit) tensors (a'b'), (a'b'c'), (a'b'c'd') and (a'b'c'd'e'), 
respectively, and similarly in a, b, c, d, e indices in Y diagonal basis. The states of the A 
and A have additional properties, 

^^aia2^3^4^5fl6^7^8^9^10 ^Gt6 ct 7 ct 8 ct 9 ct 10 Aa 1 a2a3a4(l5? (-^■•^) 

that allow one to project out the A and A states, respectively, from the 256 antisymmetric 
states (abcde). The explicit expressions for antisymmetric tensors are, for example, 

(ab) = ab — ba, 

(abc) = abc + cab + bca — bac — acb — cba (A. 9) 

etc. Important relations are 



(12) = 


-i(l'2' 


(34) = 


-i(3'4' 


(56) = 


-z(5 / 6' 


(78) = 


-z(7 / 8 / 


(90) = 


-i(9'0' 



(A.10) 
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Symmetric E = 54 dimensional space is spanned by traceless symmetric states {a'b'} = 
a'b' + b'a! (a', b' = V, 2', • • ■ , 9', 0') and J2 a ' c "' { a ' a '} with S a ' c ^ = °- Also ' important 
relations are 

{12} = l'l' + 2'2', 
{34} = 3'3' + 4'4', 
{56} = 5'5' + 6'6', 
{78} = 7V + 8'8', 

{90} = 9'9' + 0'0'. (A.ll) 

SO(10) invariants are denned in the fundamental SO (10) basis 1', 2', • • • , 9', 0' and in the Y 
diagonal basis 1, 2, • • • , 9, 0. We give some typical examples: 

$ 2 = ^a'b'c'd'^a'b'c'd' = §abcd$abcd 

$ 3 EE $ a wA'&'e'/'$cW/' = ^a&Ae/^cde/* ( A - 12 ) 
$AA EE $ a 'b'c'd>A a > b >e'f>g'A c >d>e'f>g> = ^abcd^abefg^ cd ejg, ^tc. 

Here a', 6', c', ■ • • (a, 6, c, ■ • • ) run over all the SO (10) vector (Y diagonal) indices. 

Now, the Higgs fields A, E, A, A, and $ contain 8 directions of singlets under the G 32 i 
subgroup (see Appendix of pD|). The corresponding vacuum expectation values (VEVs) are 
defined by: 

2 

(A) = ^AiA, (A.13) 

(E) = EE, (A.14) 
(A) = v r vr, (A.15) 
(A) = UrUr, (A.16) 

3 

($) = (A.17) 

i=i 
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where unit directions Ai (i — 1,2), E, vr, vr and (i — 1, 2, 3) in the F diagonal basis are: 





_ 7(1,1,0) _ 
- A (l,l,3) ~ 


J f19 4- 


(A 18"! 


Ao 


_ 4(1.1.0) _ 

^(15,1,1) 


— (56 + 78 + 90) 


(A 19) 


/? 

11/ 


9(1,1,0) 
~ ^(1.1.1) ~~ 


^ ro^fio 1 OA] t)v[« 1 7C 1 nnll 

y60 


(A 901 


vr 


_ a(i.i.O) _ 

(10,1,3) 


W >' 


( A 01 ^ 
(A.Z1J 




-(1.1,0) 






vr 


= A (10,l,3) = 


VIlo' 13579 '- 


(A.22) 


$i 


^(1,1,1) 


-V24 (1234) ' 


(A.23) 


$2 


^(15,1,1) 


^(5678 + 5690 + 7890), 

v 72 


(A.24) 


$ 3 


^(15,1,3) 


--^-([12 + 34][56 + 78 + 90]). 


(A.25) 



The upper and the lower indices indicate the SU(3) C x S£/(2) L x C/(l)y, 517(4) x SU(2) L x 
SU(2) R quantum numbers, respectively in the case of double indices. A word about notation: 
the square brackets are used for grouping of indices. This grouping of indices is used to 
emphasize the SU(2)l and SU(3)c structures within the state vectors. The square brackets 
satisfy usual distributive law with respect to summation of indices and tensor product of 
indices, e.g. 

([12 + 34][56 + 78 + 90]) 

= (1256 + 1278 + 1290 + 3456 + 3478 + 3490) 
= (1256) + (1278) + (1290) + (3456) + (3478) + (3490), 
([1, 3][5[78 + 90], 7[56 + 90], 9[56 + 78]]) 
= (1578 + 1590, 1756 + 1790, 1956 + 1978, 
3578 + 3590, 3756 + 3790, 3956 + 3978) 
= (1578, 1756, 1956, 3578, 3756, 3956) 

+ (1590,1790,1978,3590,3790,3978). (A.26) 
The unit directions appearing in VEVs satisfy the following orthonormality relations 

(*,J = 1,2), 

^ = 0, 

(i,j = 1,2,3). (A.27) 
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A- 


■A, 


= Sij 




E 2 


= 1^ 


Vr 


■ vr 


= vr 


vr 


■vr 


= 1, 




.<&.- 


= 



Due to the D-flatness condition the absolute values of the VEVs, v r and vr are equal, 



\vr\ = \vr\. 



(A.28) 



More details should be reffered to the original paper [60 
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